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Preface 


The  purpose  of  this  thesis  was  to  develop  and  validate  a  computer 
model  of  strip  geometry  explosive  flux  compression  generator.  The  Air 
Force  Weapons  Laboratory  had  designed  a  strip  generator  and  conducted 
preliminary  familiarisation  experiments  with  it.  A  computer  program  to 
predict  generator  output  was  developed.  This  computer  model  was  then 
validated  using  the  AFWL  strip  generator. 

We  would  like  to  thank  our  advisor,  Major  Timothy  L.  Skvarenina 
of  the  Air  Force  Institute  of  Technology  at  Wright-Patterson  AFB,  Ohio 
for  his  assistance  during  this  project.  We  are  especially  grateful  to 
Dr.  Robert  E.  Reinovsky  of  the  Air  Force  Weapons  Laboratory  at 
Kirtland  AFB,  New  Mexico.  His  guidance  and  assistance  were  extremely 
helpful.  We  are  also  grateful  to  Wally  Kaiser  of  Maxwell  Labs  at 
Kirtland  AFB,  New  Mexico.  His  technical  knowledge  and  enthusiastic 
support  made  the  experiments  both  enjoyadle  and  highly  productive. 

This  effort  was  jointly  funded  by  the  Air  Force  Weapons 
Laboratory  and  the  Defense  Nuclear  Agency. 
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Abstract 


A  computer  program  to  pradicc  the  performance  of  a  atrip  geometry 
explosive  magnetic  flux  compression  generator  (MCG)  ves  developed  and 
validatad. 

The  program  contains  five  separate  models  of  the  MCG,  two 
theoretical  and  three  empirically  based.  The  first  model  is  a 
theoretical  lossless  flux  conserved  model.  The  second  model  accounts 
for  resistive  losses  due  to  plate  resistance  and  a  flux  diffusion 
equivalent  resistance.  The  first  empirical  model  aodela  the  remaining 
losses  as  a  slipping  loss.  The  second  empirical  model  adjusts  the 
flux  conserved  model  with  an  additional  time  varying  inductance.  The 
final  empirical  model  is  a  simple  flux  trapping  efficiency  model. 

Actual  MCG  experiments  were  conducted  u  ing  the  Air  Torce  Weapons 
Laboratory's  strip  MCG,  The  data  from  these  experiments  was  used  to 
validate  the  computer  program.  Input  currents  varied  from  294.0  KA  to 
1100.0  KA  with  resulting  output  currents  of  2210.0  KA  to  10261.0  KA 
and  current  amplification  factors  between  6.29  and  15.75. 


COMPUTER  MODELING  AND  EXPERIMENTAL  VALIDATION 
OF  LOSSES  IN  A  STRIP  GEOMETRY  EXPLOSIVE 
MAGNETIC  FLUX  COMPRESSION  GENERATOR 


I.  Introduction 

Magnetic  flux  compression  generators  (MCGs)  are  devices  which 
create  ultrahigh  magnetic  fields  by  compressing  magnetic  fields 
between  two  metal  surfaces  to  achieve  field  multiplication.  This  field 
multiplication  is  achieved  by  converting  chemical  energy  into 
electrical  energy. 

Background 

MCGs  have  been  under  investigation  since  the  early  1950s  (1:10). 
Initial  experiments  were  done  in  strip  geometries,  explosively  driving 
two  plates  together.  Because  experiments  in  coaxial  geometry  proved 
more  efficient,  attention  soon  focused  on  coaxial  geometry  (2:5191). 

For  many  years,  the  focus  remained  on  che  cc axial  flux  compression 
generators. 

Recently,  attention  has  returned  to  the  strip  geometry  MCG. 

Because  of  their  higher  current  densities  and  relatively  simple 
construction,  the  military  is  now  considering  strip  geometry  MCGs  as 
a  possible  power  supply  for  some  of  its  advanced  weapons  and 
experisients  (3). 

By  their  nature,  explosive  MCGs  are  one  shot  devices.  For  this 
reason,  computer  models  predicting  the  MCG's  output  characteristics 
are  beneficial  in  terms  of  cost  and  time  of  development. 


Problem 


The  purpose  of  this  thesis  wes  to  develop  end  velidete  e  simple 
computer  model  of  e  strip  geometry  MCG's  operation.  The  computer  model 
accounts  for  the  major  physical  processes  involved  in  the  MCG's 
operation.  Model  validation  vas  done  using  the  Air  Force  Weapons 
Laboratory's  (AFWL)  MCG  shown  i  Figure  1-1. 


Figure  1-1.  Air  Force  Weapons  Laboratory's  MCG 


Scope 

The  study  done  for  this  thesis  was  limited  to  the  design  and 
validation  of  a  computer  model  of  a  strip  geometry  MCG.  The  loss 
mechanisms  considered  were  flux  diffusion,  plate  resistance,  sliding 
contact,  and  incomplete  compression. 

The  criteria  used  to  evaluate  the  model  were  the  ability  to 
predict  the  output  current  as  a  function  ;  time  and  the  time  to  peak 
current  amplification. 


Assumptions 

In  order  to  model  the  MCG's  operation  several  basic  assumptions 
were  made.  These  assumptions  include: 

1)  The  electrical  circuit  used  in  the  experiments  can  be  modeled 
by  the  lumped  parameter  elements  of  inductance,  capacitance,  and 
resistance  (4:327);  (This  allows  the  MCG's  operation  to  be  described 
by  RLC  circuit  analysis.); 

2)  The  magnetic  fields  are  spatially  uniform  within  the  MCG ; 

3)  The  load  and  throat  inductance  of  the  MCG  are  constant 
with  time; 

4)  The  metal  plates  are  incompressible  (5:364); 

5)  The  explosive  process  brings  the  two  plates  together  in  a 
uniform  notion  so  that  the  compression  can  be  considered  linear  (6:7); 

6)  During  the  compression  process,  the  metal  plates  offer  no 
resi'jance  to  plate  motion  (6:3);  and 

7)  Specific  heat  remains  constant  (5:362). 
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General  Approach 

Computer  Model.  In  order  to  model  the  generator's  operation  with 
an  equivalent  electrical  circuit,  basic  relationships  which  describe 
the  physical  processes  involved  had  to  be  developed.  The  first  problem 
was  to  determine  how  the  explosive  energy  was  transferred  to  the  metal 
plates.  The  detonation  of  the  explosive  accelerated  the  metal  plate 
converting  the  chemical  energy  of  the  explosive  to  kinetic  energy  or 
changing  velocity.  Although  complex  computer  codes  exist  to  describe 
this  energy  conversion,  the  Gurney  method  provides  a  quick,  simple 
approximation  for  the  changing  metal  velocity  (6:3). 

The  second  physical  process  to  be  accounted  for  was  the  time 
varying  inductance.  The  cnanging  metal  velocity  predicts  the  motion 
of  the  metal  plates  which,  in  turn,  is  used  to  calculate  the  changing 
area  enclosed  by  the  MCG.  From  this  changing  area,  a  time  varying 
inductance  was  calculated.  Combined  with  the  lumped  element  parameters 
determined  from  the  experimental  circuit,  a  totally  flux  conserved 
(lossless)  model  was  developed.  As  the  plates  are  driven  together, 
the  inductance  decreases,  converting  kinetic  energy  to  electric  energy. 

Using  this  flux  conserved  model,  the  flux  diffusion  losses  were 
calculated.  These  losses  were  modeled  by  a  time  varying  effective 
resistance.  Losses  due  to  the  copper's  resistance  were  also  computed. 
This  resistance,  combined  with  the  effective  resistance  which  accounts 
for  flux  diffusion  losses,  was  used  to  produce  a  total  resistance 
model  of  the  MCG. 

During  experimental  validation,  it  was  discovered  that  other 
physical  loss  mechanisms  were  present.  These  losses  had  to  be 
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accounted  for  in  order  to  completely  describe  the  MCG's  operation. 

Three  empirical  models  were  developed.  The  first  model  which  accounts 
for  slipping  losses  was  developed  in  a  manner  similar  to  that  of 
Knoepfel  (1:189-191).  A  second  model  which  allows  for  compression  at 
rates  faster  or  slower  than  predicted  by  the  Gurney  velocity  was  also 
developed.  A  final  model  which  is  based  totally  on  empirical 
measurements  accounts  for  a  flux  compression  efficiency. 

A  detailed  description  of  the  computer  models  can  be  found  in 
Chapter  III. 

Experimental  Validation.  Validation  of  the  computer  program 
consisted  of  performing  actual  strip  MCG  experiments.  The  experimental 
circuit  was  set  up  as  shown  in  Figure  1-2.  Constant  ircuit  parameters 
such  as  line  inductance  and  resistance  were  determined  by  discharging 
the  capacitor  bank  through  a  shorted  line  as  well  as  through  the 
AFWL  MCG.  Current  measurements  were  used  to  determine  these  values. 

Several  parameters  were  varied  to  check  the  validity  of  the 
computer  model.  These  included  changing  the  load  impedance,  the 
tamping,  the  amount  of  explosive,  a..d  the  input  voltage  to  the 
capacitor  bank. 

Five  MCG  experiments  were  successfully  performed.  Test 
measurements  included  voltages  across  the  switch  and  load,  and  current 
derivatives  through  the  switch,  line,  and  load.  Currents  were 
obtained  by  actively  integrating  the  current  derivatives.  Data  was 
displayed  on  oscilloscopes  photographed  with  polaroid  cameras. 


5 


Figure  1-2.  Experimental  Set  Up 
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A  detailed  description  or  the  experimental  apparatus  can  be  found 
in  Chapter  IV.  Chapter  V  describes  the  experimental  procedures  used. 
Finally,  the  results  comparing  the  experimental  data  with  the  computer 
models,  conclusions,  and  recommendations  can  bs  found  in  Chapter  VI. 
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determined  by  the  particular  uae  of  the  generator.  The  models  all 
describe  the  generator  inductance  as  a  time  varying  inductance.  They 
differ  in  their  description  of  the  resistance;  some  describe  it  as  a 
time  varying  resistance  while  others  describe  it  as  a  constant. 

Time  Varying  Inductance .  Inductance  varies  as  a  function  of  the 
dimensions  of  the  system.  During  the  compression  process,  the  metal 
plates  are  explosively  driven  together.  This  drastically  alters  the 
dimensions  of  the  system  and  hence  the  system  inductance.  Thrcugh  a 
change  of  variables,  this  dimensionally  changing  inductance  can  be 
modeled  as  a  time  varying  inductance  (7:147;8:334) . 

Time  Varying  Resistance.  Experts  generally  agree  that,  as  a 
minimum,  resistance  is  a  function  of  resistivity,  temperature,  and 
path  length  (8:334 ;9:S54).  Resistivity  and  temperature  are 
interrelated.  As  long  as  the  metal  plates  do  not  liquify,  resistivity 
can  be  described  as  a  linear  function  of  temperature  (10:393). 

p  -  Pn(l  +  k0)  (2-1) 

where 

Po  is  the  resistivity  of  the  plate  at  room  temperature 
k  is  a  constant  equal  to  10  /°C 

0  is  the  temperature  in  °C 

A  method  must  then  be  determined  to  describe  how  temperature  changes 
with  a  known  variable.  Because  the  thermal  diffusion  process  is  much 
slower  than  the  magnetic  diffusion  process,  temperature  can  be 
described  strictly  as  a  function  of  the  magnetic  flux  density  (1:74). 


9 


y  ■'nuimnn^ * 

a 


r.WVra’A’W  VS  V  v  v 


r.  VV^"** V ’’?r " 7  7  T  ~*r, 


9  -  x®2 


(2-2) 


where 


X  is  a  constant  equal  to  1170°C/(KG)2 
B  is  the  magnetic  £lux  density  in  megagauss  (10:392) 

As  MCG  operation  progresses,  the  magnetic  flux  is  compressed.  (A 
detailed  discussion  of  how  the  magnetic  flux  is  compressed  is  contained 
i  Chapter  III.)  This  means  that  the  magnetic  flux  density  and  hence 
th^  resistivity  can  be  made  a  function  of  time  through  a  change  of 
variables.  Additionally,  the  path  through  which  the  current  flows 
changes  during  MCG  operation.  By  using  the  familiar  expression 


(2-3) 


where 


p  is  resistivity 
£  is  path  length 
A  is  cross  sectional  area 

The  resistance  of  the  plates,  Rp,  can  be  defined  as  a  function  of  time. 

The  process  of  magnetic  flux  diffusion  is  modeled  in  the 
equivalent  electrical  circuit  approach  as  a  time  varying  resistance 
(1:185).  By  using  a  shin  layer  method,  a  time  varying  resistance 
modeling  flux  diffusion  is  defined  as 


1 

an/  T 

OoU 


(2-4) 


where 

Oo  is  conductivity 

y0  is  magnetic  permeability 

T  is  a  characteristic  time  defined  by 

1  "  (2“5) 
dt 

Because  t  is  a  function  of  time,  magnetic  flux  diffusion  can  be 
modeled  as  a  time  varying  resistance. 

Constant  Resistance  Model.  Experts  who  model  the  resistance  as  a 
constant  during  MCG  operation  roalise  that  resistance  does  in  fact 
change  with  time.  These  changes  ar<;  small  and  linear.  Thus,  they  can 
be  modeled  as  an  average  resistance  (1:83) ;ll:DI-8-3) .  In  this 
manner,  resistance  can  be  modeled  as  a  constant. 


Empirical  Modeling 

Empirical  models  are  based  on  experimental  data.  Knoepfel 
suggests  that  a  loss  mechanism  not  modeled  by  previous  methods  is  a 
slipping  contact  loss.  Slipping  losses  occur  when  plate  compression 
is  incomplete  and  flux  is  trapped  between  the  plates.  For  the  AFWL 
generator,  this  slipping  loss  occurs  on  the  surfaces  shown  as  a  drk 
line  in  Figure  2-2. 


Figure  2-2.  Slipping  Loss  Surfaces 
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Along  these  surfaces,  ’.»e  empirical  model  claims  flux  is  lost  at  a 
rate  of 

-  -2sv^ioH  (2-6) 

where 

s  is  a  mean  depth  of  ripples  in  the  slipping  surfaces 
is  t.  e  detonation  velocii  ’  shown  in  Figure  2-2 
Mo  is  the  magnetic  permeability 
H  Ls  the  magnetic  field 

In  a  manner  similar  to  that  of  Knoepfel,  it  has  been  determined  tnat 

^--(CU))1  d  (2-7) 

lo 

where 

lo  is  the  initial  current  in  the  generator 
^(t)  is  the  flux  comp re  iion  factor  determined  in  Chapter  III 
a  is  the  mean  depth  of  the  ripples 

d  is  the  initial  plate  separation 

In  this  manner,  an  output  current  can  be  determined. 

A  second  empirical  method  to  account  for  additional  losses  adjusts 
the  time  'trying  inductance  of  the  generator  predicted  by  the 
resistance  model.  This  additional  inductance  Lex  is  added  to  the  total 
time  varying  generator  inductance. 

^Gen^*)  "  1*g(0  +  Lex(t)  (2—8) 

where  I>Q(t)  is  the  time  varying  inductance  predicted  by  the  Gurney 
method.  Lex(t)  can  be  either  negative  or  positive  in  value.  A 
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negative  value  indicate*  a  faster  flux  compression  and  smaller  time 
varying  inductance.  A  positive  value  indicates  a  slower  flux 
compression  and  larger  time  varying  inductance. 

Finally,  the  generator's  performance  can  be  empirically  modeled 
by  calculating  a  flux  trapping  efficiency. 


•  Meal  Flux  Present  -  Flux  Lost  X  10GX 


Ideal  Flux Present 


(2-9) 


The  ideal  flux  present  is  the  initial  flux  present  when  generator 
operation  begins.  The  flux  lost  is  all  flux  lost  regardless  of  the 
loss  mechanism. 
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III.  Computer  Model 


Five  baeic  models  were  developed  to  described  ICG  operation.  Two 
of  these  models ,  the  flux  conserved  model  end  the  time  varying 
resistance  model  are  theoretical  equivalent  circuit  models.  The 
remaining  models,  the  sliding  contact  loss  model,  the  changing  area 
lose  model,  and  the  efficient  model  are  empirical  models. 

All  models  are  derived  from  the  flux  conserved  model.  The  models 
are  capable  of  describing  operation  for  both  rectangular  and 
trapesoidal  geometry  atrip  generators.  A  detailed  description  of  the 
program  inputs,  outputs,  and  the  program  itself  can  be  found  in 
Appendix  A. 

Theoretical  Models 

Flux  Conserved  Model.  The  flux  conserved  model  uses  an  equivalent 
electrical  circuit  approach  to  model  the  MCG's  operation.  This  model 
assumes  that  there  are  no  flux  diffusion  losses  and  that  the  metal 
plates  are  perfectly  conducting  so  that  the  circuit  has  essentially  no 
resistance.  The  circuit  relationship  is  then 

j£<LI)  -  0  (3-1) 

The  current  solution  at  any  time,  t,  is  then  given  by 

I(t)  •  (3-2) 
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L(0)  and  Io  are  initial  conditions  at  t«0 

L(t)  la  a  time  varying  ays tan  inductance 

With  no  losses,  the  MCG  operation  may  be  modeled  completely  by  a 
time  varying  inductance.  For  this  particular  circuit,  there  are  three 
inductances,  the  generator  inductance,  Lq,  the  load  inductance,  L^, 
and  the  throat  inductance,  &o .  The  throat  inductance  is  the  inductance 
of  the  small  transmission  line  connecting  the  generator  to  the  load. 
Both  and  &0  are  assumed  to  be  constant. 

The  generator  inductance  varies  with  time.  The  change  in 
inductance  results  from  the  conversion  and  partial  transfer  of  chemical 
energy  of  the  explosive  to  the  kinetic  energy  of  the  top  metal  plate. 
Rather  than  using  a  complex  computer  code  to  derive  this  energy 
transfer,  this  computer  model  uses  the  Gurney  method  vhich  provides  a 
quick,  accurate  approximation  of  the  energy  transfer. 

The  Gurney  method  is  based  on  energy  and  momentum  balances  and 
can  be  used  to  describe  a  large  variety  of  gemoetries.  For  the  strip 
generators,  the  top  plate,  explosive,  and  tamper  can  be  modeled  as  an 
asymmetric  sandwich  as  shown  in  Figure  3-1. 


Figure  3-1.  Asymmetric  Sandwich 
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In  this  particular  configuration,  the  velocity  imparted  to  the 
is  described  as 


>tei 


V  -  r/TS 
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vhere 

/TE  is  the  Gurney  energy  for  th*  particular  explosive 
—  is  the  tamper  mass  to  explosive  mass  ratio 
M 

£-  is  the  metal  mass  to  explosive  mass  ratio 

1  M 

A  is  a  constant  defined  by  *  £  (6:8) 

: — ?r 
1  *  c 

As  the  explosive  burns,  the  top  plate  is  driven  to  the  bottom  plate 
as  shown  in  Figure  3-2. 


The  angle.  A,  shown  in  Figure  3-2  is  known  as  the  Gurney  angle.  The 
Gurney  angle  is  defined  by  the  following  equation 


(3-4) 


where 

v  is  the  metal  velocity  from  Eq  (3-3) 
m 

is  the  detonation  velocity  of  the  explosive  (6:27) 

Using  the  detonation  velocity,  the  Gurney  angle,  and  simple 
geometry,  a  time  varying  area  can  easily  be  calculated.  Generation 
operation  ceases  when  the  top  plate  of  the  generator  touches  the  jet 
catcher  as  shown  in  Figure  3-3. 


Figure  3-3.  Generator  Cessation 


Using  this  time  varying  area,  the  time  varying  inductance  is  easily 
calculated  using 


L(t)  -  1  'WWk  pH  (3-5 

where 

2 

A(t)  is  the  time  varying  area  in  m 
W  is  the  plate  width  in  m 

Because  the  program  allows  for  the  width  of  the  plate  to  vary,  L(t)  is 
calculated  incrementally. 

Once  the  time  varying  generator  inductance  is  calculated,  the 
system  inductance  can  be  defined  as 


L(t)  -  Lg(t)  +  Jlo  +  1^  (3-6) 

The  flux  conserved  solution  to  Eq  (3-2)  is  then 

L_(0)  +  &o  +  K 

IIdeal<t)  "  Io  t'G CtT  VToT  1^  (3“7) 

From  this  ideal  current  the  magnetic  field  and  magnetic  flux  density 
are  easily  found  using 


H(t) 


IIdeal(t) 

W 


(3-8) 


and 


B(t)  -  UoH(t) 


(3-9) 


Time  Varying  Resistance  Model.  The  time  varying  resistance  model 
also  uses  an  equivalent  electrical  circuit  approach.  It  assumes  the 
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only  resistive  mechanism*  are  flux  diffusion  and  plate  resistivity 
The  circuit  relationship  is  then 


^<LI)  +  IR  -  0 


(3-10) 


The  current  solution  at  any  tine  is  then 


Kt) 


(3-11) 


where 

L(0)  and  Io  ere  initial  conditions  at  t-0 
L(t)  is  the  time  varying  system  inductance 
R(t)  is  the  time  varying  system  resistance 
L(t)  has  been  described  in  the  flux  conserved  model  and  the 
mechanism  for  its  generation  is  identical  in  the  time  varying 
resistance  model.  In  this  particular  model ,  it  is  assumed  that  all 
resistive  mechanisms  are  described  by  the  resistivity  of  the  copper 
plates  or  as  an  equivalent  resistance  due  to  flux  diffusion. 

Plate  resistance  is  described  by  the  well-known  relationship 

R  -  (3-12) 

p  A 

where 

p  is  resistivity 
l  is  path  length 

A  is  cross  sectional  area  in  which  the  current  flows 
As  metal  heats,  its  resistivity  changes.  As  long  as  the  plates 
remain  solid,  this  change  can  be  approximated  by  the  linear  function 


p  -  p0(l  +  ke) 
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(3-13) 


where 


Po  i ■  resistivity  at  room  temperature 
k  is  a  proportionality  constant  (~10"/°C) 

8  is  the  temperature  in  °C  (10:393) 

A  method  must  then  be  used  for  describing  how  the  temperature  changes 
with  time.  A  reasonable  approximation  for  temperature  below  melting 
point  is 

6  -  X*2  (3-14) 


where 

X  -  1170°C/(MG)2 

B  is  the  magnetic  flux  density  in  megagauss  (10:392) 

X  is  a  function  of  the  apecific  heat,  Cp,  and  the  magnetic 
permeability,  p0*  Resistivity  can  then  be  described  as 

p(t)  -  po(l  ♦  kxB2)  (3-15) 

The  path  length  also  varies  with  tine.  As  the  generator  is  compressed, 
the  path  length  decreases.  By  using  simple  geometry  and  the  Gurney 
angle,  a  time  varying  path  length  is  calculated.  The  cross  sectional 
area  of  the  plate  is  defined  as 


A  -  Wd 


(3-16) 


where 

W  ia  the  width  of  the  generator 

d  is  the  depth  of  penetration  of  the  magnetic  field 
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d  it  defined  ee 


d  -  2/cJt  (3-17) 

where 

Ko  ■  j”  end  p  is  defined  in  Eq  (3-13) 

C  is  tine 

Coupling  ell  these  tine  varying  functions  together  rest  -a  in  a  piece 
wise  continuous  tine  varying  Rp(t). 

During  generator  operation,  a  portion  cf  the  flux  diffuses  into 
the  conductor.  By  using  the  skin  layer  method  described  in  Knoepfel 
(1:184-186),  the  flux  resistance  could  be  modeled  as  a  time  varying 
resistance  by 

8^  ■  (3-18) 

*♦ 

where 

p  is  the  resistivity  described  by  Eq  (3-13) 
s^  is  a  magnetic  skin  depth 
s^  is  defined  as 

■  /k'oT  (3-19) 

where 


D 

t  *  <Hi"  ^ased  on  ***«  ideal  magnetic  fields 
It 


Coupling  these  tin*  varying  functions  together  results  in  a  time 
varying  Rd(t).  The  total  system  time  varying  resistance  is  then 
described  by 


R(t)  -  R  (t)  +  R.(t) 
P  « 


(3-20) 


The  resistance  model  solution  to  Eq  (3-10)  is  then 


(L  (0)  ♦  11.  ♦  1,1  r  /•  R  (t)  ♦  R  (t) 

ynr  *  4,  j1'  **>[/  tjt)  7T7  ♦  ^  Jtj  <3*3 


Integration  is  done  using  the  trapesoidal  approximation.  The  magnetic 
fields  are  found  in  a  manner  identical  to  that  of  the  ideal  case 
•ub.ti tutin«  IRea  tor  IId..r 

Empirical  Models 

Empirical  models  were  developed  when  it  was  discovered  that  the 
resistive  loss  mechanisms  were  not  the  primary  loss  mechanism.  Each 
of  the  following  models  uses  a  different  approach  to  attempt  to  model 
the  additional  loss  mechanisms  occurring  in  the  generator. 

Sliding  Contact  Model.  In  a  manner  similar  to  that  of  Knoepfel 
(1:189-191)*  an  empirical  model  for  slipping  contact  losses  was 
developed.  Knoepfel  accounts  for  all  losses  with  his  sliding  contact 
loss.  The  sliding  contact  model  developed  herein  assumes  the  loss 
smchanisma  are  resistive  losses  and  sliding  contact  losses.  The 
sliding  contact  loss  mechanism  is  due  to  flux  trapping  when  two  plates 
are  driven  together  and  the  plate  surfaces  do  not  come  together 
smoothly. 


It  was  assumed  that  tha  sliding  contact  current  could  be 


described  by 


ISlip(t) 


W*> 


1  - 


2s 


(3-22) 


where 


Io 


is  the  initial  current 


IRcs(t)  is  the  tine  varying  current  using  tha  resistive  model 
s  is  a  mean  ripple  depth 

d  is  the  initial  plate  separation 

Solving  Eq  (3-21)  for  a  yields 


d 

*  "  I 


(t).  -  log  ♦  i 


) 


(3-23) 


Usirg  experimental  data,  s(t)  was  calculated  for  each  test  using 
the  actual  current  in  place  of  Islip(t) .  These  curves  were  linearised 
using  three  tines  during  generator  operation.  These  times  are  shown 
in  Figure  3-4.  The  first  natch  point  corresponds  to  when  the  top  plate 
has  completed  its  run  dov j  the  cutting  wedge  and  begins  to  compress 
against  the  bottom  plate.  The  second  match  point  corresponds  to  when 
the  top  plate  of  the  generator  is  about  to  begin  its  run  down  the  end 
wedge.  The  final  match  point  occurs  when  the  maximum  experimental 
output  current  is  reached. 

The  ripple  depth  at  each  match  point  is  calculated  from  Eq  (3-23), 
and  the  ripple  depth  at  any  point  between  the  match  points  is  found 
through  linear  interpolation.  From  the  sets  of  data,  best  and  worst 
case  slipping  losses  were  modeled.  Igiip^  *•*  ':ben  calculated  using 
Eq  (3-21)  and  the  experimentally  derived,  computer  linearised  ripple 
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depth,  s(t).  The  magnetic  field*  are  again  found  in  a  manner  identical 
to  that  of  the  ideal  case  substituting  IIdeal  with  1^.^. 


Point  1 


r'X 


Point  3 


Figure  3-4.  Slipping  Model  Match  Points 


Changing  Area  Loss  Model.  Since  the  ripple  depth,  s,  could  not  he 
tied  to  any  known  parameter,  a  changing  area  loss  model  vas  developed. 
This  model  assumes  that  the  changing  inductance  derived  from  the 
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Gurney  approximations  may  ba  incorrect.  It  assumes  that  the  generator 
inductance  can  be  defined  as 


W!>  ■  V*>  ♦  L„<t) 


(3-24) 


where 

Lg(t)  is  the  generator  inductance  calculated  from  the  Gurney 
approximations 

LQX(t)  is  an  additional  inductance  caused  by  incomplete  Gurney 
compression 

If  we  replace  Lg(t)  in  Bq  (3-21)  with  L-  (t),  the  principal  effect 
is  to  decrease  the  flux  compression  factor,  y(t). 


Lr(0)  +  f>o  ♦  L. 

T<t>  ’ 


(3-25) 


Because  Lg(t)  is  determined  by  a  changing  area,  so  too  is  (t)» 
This  model  assumes  that  the  changing  area  current  can  be  described  by 


^area*^ 


l*g(0)  +  to  + 

10  LG(t)  ♦  L#xrty  ♦  I#  ♦  ^ 


A')  ♦  vtJ 


dt 


(3-26) 


Solving  Bq  (3-26)  for  L  (t)  yields 


L  (0)  +  A0  +  L. 

W*>  -  *• 


Cares 


TtT 


exp 


|‘/r 


dt  -  (I*6(t)  + 


u  *  h,} 


(3-27) 
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Using  experimental  data,  L^(t)  was  calculated  for  each  taat  using 
tha  actual  currant  in  placa  of  ICarca(t)*  °QC*  calculated,  a  changing 
area  was  calculated  using  Eq  (3-5). 


L  (t>« 

ex 


(3-28) 


This  changing  area  was  than  graphed  for  analysis. 

Efficiency  Modal.  In  a  manner  similar  to  that  of  Fowler,  a  flux 
trapping  efficiency  model  was  developed  (12:26).  Using  the  well-known 
relationship  for  flux  of 


$  •  LI 


(3-29) 


a  time  varying  $  was  developed.  In  this  model,  it  is  assumed  that 
L„(t)  from  the  Gurney  calculations  is  a  valid  representation  of  the 
time  varying  inductanca.  This  model  makes  no  attempt  to  physically 
explain  where  or  why  the  flux  loss  occurs.  It  models  all  losses, 
whether  due  to  slipping  or  due  to  inaccuracies  in  the  Gurney  derived 
inductance  aa  a  time  varying  flux  trapping  efficiency.  This  flux 
trapping  efficiency  is  defined  by 


(3-30) 


where 

4>o  is  the  initial  flux  in  the  system 
$(t)  is  calculated  from  test  data 
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This  efficiency  it  smoothed  out  end  linearised  and  a  beat  and  worst 
cate  effieianeiaa  are  determined.  An  efficiency  current  ia  then 
determined  using 

W'>  ‘  ClR..(t)  <*-»> 

The  magnetic  fields  are  found  in  a  manner  identical  to  that  of  the 
ideal  case  substituting  with  . 
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Experiments  to  validate  the  computer  model  mere  conducted  et  the 
Air  Force  Weapon*  Laboratory  explosive  test  facility  located  on  the 
McCormick  Range  at  Kirtland  AFB,  Hem  Mexico.  The  major  components  and 
diagnostic  equipment  used  in  the  magnetocuaulative  generator  tests  are 
described  in  detail  in  this  chapter.  The  major  components  consist  of 
the  capacitor  bank,  the  solid  dielectric  explosive  closing  switch, 
the  fire  control  system,  the  transmission  line,  and  the  MCG  and  load 
coil  (See  Figure  4-1).  The  primary  diagnostic  devices  used  in  the 
experiments  mere  Rogowski  coils,  voltage  probes,  and  oscilloscopes. 


Major  Components 


Capacitor  Bank.  The  capacitor  bank  provided  the  initial  feed 
current  to  the  MCG.  The  bank  consisted  of  fifty  6.1  yF  60  KV 
capacitors  connected  in  parallel  to  yield  a  total  nominal  capacitance 
of  305  uF.  Charge  voltages  between  15  KV  and  45  KV  mere  used.  The 
bank  mas  charged  by  a  Maxwell  Laboratory  chargiug  system  and  discharged 
into  the  transmission  line  and  MCG  through  a  solid  dielectric  closing 


smitch. 


Closing  Smitch.  The  closing  smitch  used  during  the  experiments 
was  a  solid  dielectric  explosive  switch.  The  switch  consisted  of  two 
6  foot  long  by  8  inch  wide  by  0.5  inch  thick  steel  plates.  The 
connection  of  the  closing  switch  to  the  transmission  line  and 
capacitor  bank  is  shown  in  Figure  4-1. 

One  of  the  plates  was  connected  via  the  solid  inner  conductors  of 


numerous  parallel  coaxial  cables  to  the  capacitor  bank.  The  other 


Figure  4-1.  MCG  Experimental  Set  Up 


plate  wee  attached  to  the  top  of  the  transmission  line  and  could  be 
pulled  back  to  facilitate  cleaning  and  rebuilding  of  the  switch. 

Sheets  of  mylar  sandwiched  between  the  two  aligned  plates  acted  as 
the  solid  dielectric  insulation  for  the  switch.  A  1  mil  thickness  of 
mylar  was  used  for  each  1  KV  of  potential  across  the  plates.  Plate 
separation  was  determined  by  the  number  of  sheets  of  mylar  required 
for  the  test  voltage  applied. 

Before  each  test,  two  RP2  PETN  explosive  shaped  charge  detonators 
were  placed  into  0.5  inch  holes  which  were  drilled  along  the  top  edge 
of  the  plate  connected  to  the  transmission  line.  Closing  of  the  switch 
occurred  when  the  two  detonators  were  fired,  rupturing  the  mylar 
insulation  and  allowing  breakdown  of  the  switch.  The  capacitor  bank 
then  discharged  its  energy  through  the  switch  into  the  transmission 
line  and  the  MCG. 

Because  of  its  destructive  nature,  the  switch  required  rebuilding 
after  each  shot.  Rebuilding  consisted  of  replacing  the  mylar  sheets 
between  the  plates,  replacing  the  detonators,  and  reconnecting  the 
detonation  lines  to  the  X  Units. 

Firing  System.  The  firing  systems  consisted  of  two  TRW  Model  64A 
Trigger  Delay  Generators  and  a  Los  Alamos  National  Laboratory  X  Unit. 
The  trigger  delay  generators  provided  a  variable  time  delayed  output 
pulse  to  fire  the  Z  Unit  and  trigger  the  other  diagnostic  equipment. 

The  X  Unit  provided  the  current  pulse  used  to  fire  the  detonators  in 
both  the  switch  bunker  and  the  explosive  on  the  MCG. 

Transmission  Line.  The  transmission  line  consisted  of  two 
parallel  plates  of  60  mil  aluminum,  40  feet  long  and  4  feet  wide. 


This  length  allowed  the  feed  current  to  be  delivered  to  the  MCG  at  a 
safe  distance  from  the  bunker  and  switch,  avoiding  damage  to  the 
switch  during  the  MCG  explosions. 

The  two  aluminum  sheets  are  separated  by  twelve  sheets  of  5  mil 
thick  6  foot  wide  mylar.  This  allows  for  a  hold  off  voltage  of  60  KV 
at  1  mil  per  kilovolt.  The  total  insulation  thickness  of  60  mils 
results  in  a  line  capacitance  of 


C  -  eA/d 
-  0.26  pF 


(4-1) 


Since  Ck£ne  «  Cfiank,  Cj^ne  was  neglected  in  all  circuit  calculations. 
The  plate  edges  were  further  insulated  against  edge  field  enhancement 
effects  by  the  placement  of  0.5  inch  diameter  tygon  tubing  between  the 
sixth  and  seventh  sheets  of  mylar  along  the  edge  of  the  entire  length 
of  the  line  (3). 

Because  the  end  of  the  transmission  line  was  connected  directly 
to  the  MCG,  it  was  destroyed  with  each  test.  To  avoid  replacing  the 
entire  transmission  line,  6  foot  long  by  4  foot  vide  end  pieces  were 
bolted  onto  the  end  of  the  line  using  clamping  bars  to  ensure  good 
electrical  contact.  Sandbags  protected  the  clamping  bars  during  the 
explosions.  Following  each  test,  the  end  pieces  and  final  sections  of 
mylar  were  replaced. 
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Four  Rogowski  coils  were  used  for  each  teat,  laced  as  shown  in 
Figure  4-2. 

Switch  Bunker 


r-  -  -  -  -  - - 


Voltage  Load 

Probe  Voltage 

Probe 


Figure  4-2.  Placement  of  Test  Diagnostic  Devices 

Coils  A  and  B  were  used  to  measure  the  load  current  derivative. 

Coil  C  was  used  to  measure  the  feed  current  derivative  until 

crowbarring  occurs  and  Coil  D,  located  in  the  switch  bunker,  measured 

the  derivative  of  the  return  current  of  the  transmission  liue. 

Both  handvound  and  etched  Rogowski  coils  were  used.  The  handwound 

9 

Rogowski  coils  provided  a  1  volt  signal  for  every  10  A/ sec  present 
while  the  etched  Rogowski  coils  provided  a  1  volt  signal  for  every  10*^ 
A/sec  present.  Coil  A  was  etched,  while  coils  B,  C,  and  D  were 


handwound 


To  enhance  coil  survivability  during  the  explosion,  all  four 
Rogowski  coils  were  placed  inside  through  tygon  tubing.  Coils  A,  B, 
and  C  were  destroyed  during  each  test  and  had  to  be  replaced  prior  to 
the  next  shot.  Coil  D  remained  undamaged  throughout  the  explosion  and 
thus  provided  a  uniform  means  to  compare  the  input  parameters  from 
test  to  test. 

The  current  derivatives  from  coils  A,  B,  and  D  were  actively 
integrated  using  Textronixs  Type  0  Operational  Amplifier  integrators. 
The  output  from  these  integrators  was  fed  into  oscilloscopes  and 
recorded  on  film. 

Voltage  Probes.  At  least  one  voltage  probe  vas  used  in  each 
experiment;  the  one  located  in  the  switch  bunker  was  used  to  measure 
the  voltage  across  the  switch.  In  the  experiments  where  a  MCG  load 
coil  of  10.31  nH  was  used,  a  second  voltage  probe  was  connected  to  the 
load  coil  to  measure  the  voltage  across  the  load  during  the  operation 
of  the  MCG.  The  voltage  probes  consisted  of  resistive  dividers  to 
measure  voltage  and  were  connected  as  shown  in  Figure  4-2.  In  both 
voltage  probes,  the  two  resistors  were  filled  with  a  copper  sulfate 
solution  with  copper  electrodes  sealing  each  end.  The  copper  sulfate 
solution  was  adjusted  to  give  a  nominal  resistance  of  20  Kft.  Coaxial 
diagnostic  cables  fed  -he  signals  to  50  R  terminators  across  the 
oscilloscope  inputs.  Therefore,  a  2.5  volt  signal  was  measured  at  the 
oscilloscope  for  every  1  KV  across  the  switch  or  load  coil.  The  outer 
cable  shields  were  shorted  together  at  the  resistors  and  at  the 

i 

oscilloscope  inputs  to  balance  any  induced  shield  currents.  The  two 
signals  from  each  probe  were  then  connected  to  a  differential  amplifier 
which  eliminated  any  common  mode  signals  (13:182). 
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Oscilloscopes.  Data  from  each  test  vaa  recorded  on  fourteen 
Textronixs  7603  oscilloscopes.  Four  oscilloscopes  recorded  signals 
genarated  from  each  load  Rogovski  coil  with  two  recording  dl/dt  and  two 
recording  the  integrated  current  waveforms.  Two  oscilloscopes  measured 
the  differential  voltage  across  the  load  when  the  10.31  nH  load  was 
used.  One  oscilloscope  was  used  to  record  the  line  dl/dt  from  coil  C. 
Another  oscilloscope  recorded  the  line  dl/dt  from  coil  D.  Two  other 
oscilloscopes  recorded  the  integrated  current  from  coil  D  and  the 
differential  voltage  across  the  switch.  The  data  from  each  of  the 
oscilloscopes  was  recorded  by  Polaroid  cameras  which  were  triggered  by 
the  delay  generators. 


V.  Experimental  Proceduraa 


Five  atrip  geometry  magneto cumulative  generator  experiments , 
referred  to  aa  MCG  shots »  were  performed  in  order  to  obtain  data  with 
which  the  computer  smdel  could  be  compared  and  validated.  Prior  to 
every  MCG  shot,  two  preliminary  experiments  were  conducted.  The  first 
of  these,  referred  to  as  a  line  short  shot  was  required  to  obtain  ele¬ 
ment  values  for  the  equivalent  circuit  of  the  capacitor  bank  and 
transmission  line.  The  second,  referred  to  as  a  MCG  short  shot,  was 
conducted  in  order  to  calibrate  the  diagnostic  Rogowski  coils  which 
would  be  used  during  the  actual  MCG  shot.  This  chapter  describes  in 
detail  the  general  procedures  that  are  common  to  all  three  of  these 
shots.  It  then  describes  the  experimental  procedures  which  are  unique 
to  the  line  short  shot,  the  MCG  short  shot,  and  the  MCG  shot. 

General  Procedures 

Prior  to  each  shot  several  experimental  procedures  common  to  all 
shots  were  performed.  These  procedures  include  diagnostic  line  tests, 

Z  Unit  and  detonation  line  tests,  and  preparation  of  the  oscilloscopes. 

Diagnostic  Line  Tests.  To  preclude  the  loss  of  experimental  data 
due  to  faulty  diagnostic  lines,  these  lines  were  tested  prior  to  each 
shot.  A  repetitive  train  of  pulses  from  the  TRW  Deley  Generators  was 
delivered  to  the  MCG  connection  end  of  the  transmission  line  via  a  test 
line.  This  test  signal  was  connected  to  each  diagnostic  line  in  turn. 
The  return  signal  was  then  monitored  on  an  oscilloscope  to  insure  that 
the  line  was  operating  reliably  prior  to  being  connected  to  a  Rogowski 
coil  or  voltage  probe. 


X  Unit  and  Detonation  Line  Test.  Prior  to  each  testv  the  X  Unite 
and  detonation  line*  were  teated  for  reliable  operation  and  timing. 
Blown  bridge  wires  were  attached  to  the  detonation  lines  to  simulate 
detonators  during  the  test.  Using  a  pencil,  a  carbon  path  was  marked 
between  the  two  terminals  where  the  bridge  wire  was  originally 
connected,  providing  a  path  for  breakdown  when  the  X  Units  fired.  The 
resulting  flash  and  distinctive  pop  accompanying  this  breakdown  gave 
clear  indication  that  the  X  Units  were  operating  correctly. 

Preparation  of  the  Oscilloscopes.  Preparation  of  the 
oscilloscopes  prior  to  each  shot  included  calibration,  setting  of 
graticules,  and  setting  of  baselines.  Once  the  scopes  were  calibrated 
the  graticules  and  baselines  were  adjusted  to  levels  which  would 
ensure  clear  data  recording  of  measured  signals  with  cameras. 

Just  prior  to  the  shot,  the  shutters  of  the  Polaroid  cameras  were 
opened  to  record  the  graticules.  The  graticules  were  then  turned  off 
to  prevent  them  from  being  recorded  again  during  the  shot.  As  the 
capacitor  bank  was  being  charged,  the  shutters  were  opened  one  final 
time  in  order  to  record  the  date  signals.  The  TRW  Delay  Generators 
were  coordinated  with  the  data  signals  from  the  various  Rogovski  coils 
and  voltage  probes. 


Line  Short  Shots 

A  total  of  four  line  short  shots  were  performed  to  determine 
circuit  element  values  for  the  equivalent  circuit  model  of  the 
capacitor  bank,  switch,  transmission  line,  and  connections.  This 
section  describes  the  line  short  shot  and  the  general  procedures  to 
obtain  these  circuit  element  values. 


For  lino  short  shots  the  ends  of  the  transmission  line  where  the 
MCG  would  later  be  connected  were  shorted  together  giving  rise  to  the 
name  line  short  shot.  The  two  load  Rogowski  coils  and  the  end  of  the 
line  current  Rogowski  coil  were  looped  around  the  shorting  bar  at  the 
end  of  the  line. 

The  capacitor  bank  was  charged  and  then  discharged  through  the 
switch  and  the  transmission  line.  This  resulted  in  an  exponentially 
damped  sinusoidal  waveform  (See  Figure  5-1). 
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Figure  5-1.  Typical  Short  Circuit  Plot 

The  only  known  values  in  these  experiments  were  the  bank 
capacitance,  and  the  charge  voltage,  Vo.  All  other  capacitances 

were  .fleeted  because  they  were  much  less  than  CB>n^.  The  period,  T, 
was  determined  directly  from  the  oscilloscope  traces. 


The  solution  to  s  tins  invariant  SLC  circuit  is  vell-knovn  Cl* 132) . 
A  least  squares  fit  vas  dona  to  determine  the  exponential  damping 
factor,  a,  shown  in  Figure  5-1.  Knowing  T,  C,  and  a,  the  circuit 
inductance  was  calculated  using 

t  -  2nur 

-  2tt/«u>0)2  -  (o)2>1/2 

where  w0  -  ^ 
yielding 

L  -  t2/C(4tt2  +  TO2) 
and  the  circuit  resistance  vas  determined  by 

R  ■  2otL  (5-3) 

In  MCG  work,  the  input  circuit  is  crowbarred  away  from  the 
generator  after  the  first  quarte;  cycle  at  Ipeafc.  A  comparison  of  the 
values  of  Ipe«k  and  quarter  cycle  time  using  the  time  invariant 
RLC  solution  and  an  approximate  LC  solution  showed  that  R  could  be 
neglected  without  significant  loss  of  accuracy. 

MCG  Short  Shots 

This  section  describes  in  detail  the  MCG  short  shot  and  the 
general  procedures  for  calculating  the  total  system  inductance  and  the 
sensitivities  of  the  Rogovski  coils.  A  minimum  of  one  MCG  short  shot 
was  performed  prior  to  each  MCG  shot.  For  a  successful  MCG  shot,  the 
entire  system  must  operate  reliably.  When  one  part  of  the  system  was 
found  to  have  a  problem,  it  was  fixed  or  replaced  and  another  MCG 


(5-1) 


(5-2) 


short  teat  was  performed.  For  thie  reason,  as  many  as  four  MCG  short 
shots  were  performed  prior  to  one  actual  MCG  shot.  These  shots  provide 
a  final  test  of  the  entire  system  to  ensure  all  components  were 
operating  correctly. 

These  shots  provided  oscilloscope  traces  of  the  current  outputs 
from  which  the  total  system  inductance  would  be  determined.  They  also 
provided  dl/dt  traces  from  each  Rogowski  coil  so  that  their  sensitivi¬ 
ties  could  be  experimentally  determined  prior  to  the  actual  MCG  shot.. 

For  the  MCG  short  shot,  the  MCG  was  bolted  to  the  end  of  the 
transmission  line  in  the  same  manner  us  for  the  actual  shot.  The  only 
differences  were  the  lack  of  explosive  on  the  generator  and  the  place¬ 
ment  of  the  load  voltage  probe.  The  two  load  coil  Rogowski  coils  were 
placed  through  the  load  coil  and  the  end  of  the  line  return  Rogowski 
coil  was  placed  around  the  bottom  return  transmission  line's  plate  at 
the  connection  to  the  MCG.  For  the  actual  MCG  shot,  these  were  moved 
to  the  load  coil. 

The  capacitor  bank  was  charged  to  15  KV  for  each  MCG  short  shot. 
The  capacitor  bank  was  then  discharged  through  the  switch,  transmission 
line,  and  the  MCG.  This  again  resulted  in  a  damped  sinusoidal  waveform 

The  oscilloscope  traces  from  the  MCG  short  shot  were  used  in  order 
to  determine  the  total  system  inductance  of  the  capacitor  bank,  the 
switch,  the  transmission  line,  connections,  and  the  MCG.  The  period, 

T,  of  the  damped  sinusoidal  waveform  was  read  directly  from  the  traces. 

With  a  known  bank  capacitance  and  neglecting  system  resistance, 
the  total  system  inductance  could  be  determined  using 

T  -  2m4TCBank 


(5-4) 


which  whan  aolvcd  for  tha  inductanca  yialda 

(5-5) 

Oacilloacopa  tracaa  of  dl/dt  from  aach  Rogowski  coils  wera  used 
to  axparimantally  datermina  tha  sensitivity  of  the  Rogowski  coils  which 
would  be  used  during  tha  succeeding  MCG  shot.  Knowing  the  capacitor 
bank  charge  voltage ,  Vo»  the  total  inductanca,  LT,  and  neglecting 
system  resistance,  tha  maximum  dl/dt  is  given  by 


■  ft 

Measuring  tha  maximum  voltage  of  tha  dl/dt  trace,  the  sensitivity 
of  tha  Rogowski  coil  is  given  by 

Sensitivity  -  /Vmax  (5-7) 

With  the  integrator  time  constant,  TI,  the  sensitivity  of  tha 
Rogowski  coil  and  the  voltage  peak  of  the  integrated  dl/dt  trace,  Vp, 
Ipeak  was  calculated  using 

IPeak  "  yp  x  Sensitivity  X  TI  (5-8) 

This  Ipeak  was  then  compared  with  the  value  predicted  by  the  time 
invariant  RLC  solution  as  a  check  of  the  integrators.  This  procedure 
was  performed  for  each  coil  to  be  used  during  the  actual  MCG  shot. 

MCG  Shots 

A  total  of  five  MCG  shots  were  performed  to  obtain  data  for 
validation  of  the  computer  model.  This  section  details  the  set  up, 
procedures,  and  input  conditions  for  each  of  these  shots. 


40 


The  explosive  used  in  ell  shots  was  DuPont  Detashset.  As 
described  earlier  in  this  chapter,  the  total  system  inductence  as  veil 
as  the  sensitivities  of  the  Rogovski  coils  were  determined  during  a 
MCG  short  shot  just  prior  to  the  actual  MCG  shot.  The  total  inductance 
varies  with  the  length  of  the  transmission  line  end  pieces  and  the 
connection  to  the  MCG.  After  the  MCG  short  shot,  the  X  Units  and 
detonation  lines  were  once  again  tested,  the  switch  rebuilt,  and  the 
oscilloscopes  prepared.  During  this  period  of  time  explosives  are 
readied  on  the  MCG  by  contractor  personnel  from  the  Civil  Engineering 
Research  Facility  (CERF) .  Each  generator  is  identical  except  for  the 
site  of  the  load  coil.  For  three  shots  the  loed  coil  had  an  inductance 
of  10.31  nH.  For  two  shots,  a  smaller  load  was  used;  one  had  an 
inductance  of  0.2213  nH  and  the  other  had  an  inductance  of  0.5701  nH. 
The  MCG  was  placed  on  a  plate  of  one  inch  steel  which  was  separated 
from  the  MCG  by  several  sheets  of  mylar.  This  steel  plate  is  placed 
on  a  bed  of  cinder  blocks.  A  smaller  40  inch  long  by  8  inch  wide  by 
3/16  inch  thick  steel  tamper  is  placed  on  top  of  the  explosive  for 
every  shot.  The  amount  of  cement  tamper  placed  on  this  steel  sheet, 
as  well  as,  the  bank  voltage  varied  from  shot  to  shot.  Shown  in 
Table  5-1  is  a  summary  of  the  input  conditions  for  the  five  MCG  shots. 


Table  5-1 


Smeary  of  Input  Conditions 


Teat 

V(KA) 

L  (nH) 

I  (KA) 

Explo¬ 

sive 

Weight 

(lbs) 

Steel 

Weight 

(lbs) 

Cernnt 

Weight 

(lbs) 

^Load 

(nH) 

1 

44.88 

499.500 

1100.00 

1.85 

4.911 

5.0 

10.31 

2 

25.00 

437.700 

294.00 

8.16 

4.911 

20.0 

10.31 

3 

25.40 

487.693 

630.00 

5.90 

4.911 

10.0 

10.31 

4 

27.70 

540.380 

652.00 

7.37 

4.911 

20.0 

0.5701 

5 

24.99 

487.690 

620.00 

6.40 

4.911 

15.0 

0.2213 

r 


VI .  Results  and  Conclusions 

This  chapter  praaanta  tha  raaulta  of  tha  computer  models  and 
compares  than  with  tha  experimental  data.  In  addition,  it  indudea 
racoons ndat ions  for  further  atudy. 

Raaulta 

Tha  goal  of  thia  atudy  waa  to  develop  a  computer  nodal  of  tha 
atrip  geometry  MCG  and  to  validate  it  through  actual  MCG  ahota.  Tha 
primary  criterion  for  evaluating  tha  coaputer  nodal  waa  tha  nodal* a 
ability  to  reproduce  the  current  curvea  produced  by  actual  MCG 
experiments.  Accurate  reproduction  of  thaae  curvea  indicates  that  the 
conputar  program  ahould  be  able  to  reliably  predict  tha  current  curvea 
for  future  reata.  Thia  auction  ahowa  the  reaulta  of  the  actual  MCG 
ahots  and  conparea  then  with  tha  results  of  tha  coeqmter  models. 

A  total  of  aix  experiments  were  conducted  using  tha  AFVL  strip 
MCG.  Five  of  these  experiments  were  successful.  The  unsuccessful 
experiment  was  caused  by  a  destructive  failure  of  one  of  the  capacitors 
in  the  capacitor  bank.  The  results  of  the  five  successful  shots  are 
summarised  in  Table  6-1. 

A  plot  of  the  actual  current  from  the  first  experiment  is  shown 
in  Figure  6-1.  Actual  currents  from  the  remaining  tests  are  contained 
in  Appendix  B.  In  four  of  the  tests,  crowbar  occurred  at  or  near  Ip^^ 
of  the  sinusoidal  input  current  where  dl/dt  is  equal  to  sero.  In  the 
second  test,  the  generator  crawbarred  late  resulting  in  an  iuput 
current  of  294.0  KA.  Input  currents  varied  fron  294.0  to  1100.0  KA 
while  output  current  peaks  varied  fron  2210.0  to  10261.0  KA. 


Table  6-1 


Summary  of  Results 


Test 

Input 

Current  (KA) 

Maximum  Output 
Current  (KA) 
Ideal/Actual 

Current  Ampli¬ 
fication  Factor 
Ideal/Actual 

1 

1100.00 

15188.6/6928.0 

13.81/6.29 

2 

294.0 

4128.7/2210.0 

14.04/7.50 

3 

630.78 

8715.3/5355.0 

13.82/8.50 

4 

652.66 

35646.0/10261.0 

54.62/15.75 

5 

620.05 

29556.4/6031.0 

47.68/9.73 

Computer  Model  Results.  Based  on  the  input  parameters  the 
computer  program  calculates  and  prints  the  currents  predicted  by  each 
of  five  models:  the  flux  conserved  model,  the  resistance  model,  the 
sliding  contact  loss  model,  the  changing  area  model,  and  the  efficiency 
model.  The  results  of  these  models  are  best  described  in  terms  of  how 
well  they  correspond  to  the  actual  experimental  data.  The  results  are 
considered  individually. 

Flux  Conserved  Model.  The  flux  conserved  model  calculates  the 
output  current  as  if  no  losses  occur.  The  current  is  a  function  of 
only  the  time  varying  inductance,  which  in  turn  is  a  function  of  the 
detonation  velocity  aud  the  Gurney  angle.  The  time  varying  inductance 
which  corresponds  to  a  Gurney  angle  of  5.6°  from  the  first  test  is 
shown  in  Figure  6.2 

The  flux  conserved  current  versus  the  test  current  for  the  first 
teet  is  shown  in  Figure  6-3;  the  curves  for  the  remaining  tests  are 
shown  in  Appendix  B.  Points  1,  2,  and  3  from  Figure  3-4,  corresponding 


0. 00  20.00  40.00  60.00  80.00  100.00  120.00  140.00 

TIME  IN  MICROSECONDS 


to  key  generator  operation  points,  are  shown  in  this  figure.  A 
comparison  of  these  two  curves  show  that  the  flux  conserved  model 
provides  a  good  approximation  of  the  actual  currents  for  the  first  SOX 
of  generation.  Because  the  flux  conserved  model  neglects  all  losses, 
these  plots  show  that  for  the  first  50X  of  operation,  the  system  is 
basically  lossless.  During  the  last  50Z  of  operation  this  model 
predicts  higher  currents  than  those  actually  present.  This  indicates 
the  presence  of  loss  mechanisms  which  must  be  accounted  for  to  obtain  a 
complete  model  of  the  generator. 

Resistance  Model.  The  resistance  model  theoretically  calculates 
the  losses  due  to  the  plate  resistance  and  flux  diffusion  equivalent 
resistance.  A  plot  of  the  time  varying  resistance  for  the  first  test 
is  shown  in  Figure  6-4;  curves  for  the  other  tests  are  shown  in 
Appendix  B.  A  plot  of  the  time  varying  resistance  model's  current 
versus  the  test  current  for  the  first  test  is  shown  in  Figure  6-3. 

These  plots  clearly  indicate  that  plate  resistance  and  flux 
diffusion  resistance  play  only  a  limited  role  as  loss  mechanisms. 
Resistance  effects  are  negligible  during  the  initial  50%  of  generator 
operation,  again  indicating  the  system  is  essentially  lossless  during 
this  time  period.  During  the  final  SOX  of  operation,  the  time  varying 
resistance  predicts  a  current  about  IX  lower  than  that  calculated  by 
the  flux  conserved  model. 

An  examination  of  the  resistance  shows  that  the  plate  resistances 
exceed  the  flux  diffusion  equivalent  resistances  by  two  orders  of 
magnitude  throughout  generator  operation.  This  indicates  that  flux 
diffusion  plays  only  a  minor  role  as  a  loss  mechanism. 
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The  resistivity  which  would  be  required  were  resistance  alone  to 
account  for  all  losses  in  the  generator's  operation  was  calculated. 

This  resistivity  was  found  to  be  at  least  two  orders  of  magnitude 
larger  than  the  resistivity  of  vaporised  copper.  Since  the  magnitude 
of  the  ideal  B  fields  was  not  sufficient  to  cause  surface  vaporisation, 
a  resistance  this  large  was  not  reasonable.  The  resistance  model's 
poor  prediction  during  the  final  50Z  of  the  generator's  operation  shows 
that  additional  loss  mechanisms  must  be  considered. 

Sliding  Contact  Model.  Unlike  the  flux  conserved  and  resistance 
siodels,  the  sliding  contact  model  is  based  on  empirical  data.  A  time 
varying  mean  ripple  depth,  s,  was  determined  for  the  first  three 
experiments  with  a  load  inductance  of  10.31  nH.  A  plot  of  the  s  curve 
for  the  first  experiment  is  shown  in  Figure  6-5,  while  the  s  curves 
for  the  second  and  third  tests  are  presented  in  Appendix  B.  Points  1, 
2,  and  3  from  Figure  3-4,  corresponding  to  key  generator  operation 
points,  are  shown  in  this  figure. 

The  bubble  height  gives  a  measure  of  how  much  flux  is  lost  due  to 
the  sliding  losses.  Since  s  could  not  be  linked  to  any  parameter,  a 
best  and  worst  case  sliding  contact  loss  model  was  developed.  Plots 
of  the  slipping  loss  model's  best  and  worst  case  currents  versus  the 
test  current  for  the  first  test  are  shown  in  Figure  6-6  with  plots  for 
later  tests  in  Appendix  B. 

A  comparison  of  these  curves  shows  that  the  ranges  described  by 
the  curves  are  reasonable.  All  test  currents  fall  between  the  best  and 
worst  case  curves  of  the  slipping  loss  model. 

Changing  Area  Loss  Model .  Instead  of  modeling  the  additional 
losses  as  the  bubble  height,  this  model  assumes  that  the  Gurney 
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Figure  6-6.  Slipping  Current*  Teet  1 


approximations  for  tha  changing  inductance  are  not  valid.  Figure  6-7 
ahova  a  plot  of  tha  additional  area  which  must  be  added  for  the  model 
current  to  match  tha  experimental  current.  Points  1,  2,  and  3  from 
Figure  3-4  are  shown  on  this  figure.  Figure  6-8  shows  tha  Guruey 
based  area,  the  additional  area,  and  the  sum  of  these  areas;  similar 
curves  for  the  second  and  third  test  are  contained  in  Appendix  B. 

Negative  areas  indicate  that  flux  compression  is  occurring  faster 
than  that  predicted  by  the  Gurney  approximations  based  time  varying 
inductance.  This  occurs  early  during  generator  operation.  Positive 
areas  indicate  that  flux  compression  is  occurring  slower  than  that 
predicted  by  the  Gurney  approximations.  This  occurs  late  during  MCG 
operation.  A  momentum  balance  during  this  time  period  shows  that  the 
momentum  of  the  metal  plate  due  to  tha  explosive  initially  is  orders 
of  magnitude  larger  than  the  momentum  of  the  metal  plate  due  to  magnetic 
pressure.  By  the  time  the  top  plate  has  started  down  the  jet  wedge, 
the  two  momentuma  are  essentially  equal  with  the  explosive  momentum 
always  slightly  larger.  This  means  the  Gurney  approximations  are 
probably  not  valid  during  the  final  siicroseconds  of  generator  operation. 

Beat  and  worst  case  changing  area  models  were  developed.  The 
best  and  worst  case  changing  area  model  currents  for  the  first  test  are 
shown  in  Figure  6-9  while  similar  curves  for  the  remaining  tests  are 
contained  in  Appendix  B.  A  comparison  of  these  curves  shows  that  the 
experimental  currents  fall  between  the  best  and  worst  case  curves  for 
the  first  and  third  tests.  The  two  curves  are  a  bit  higher  than  the 
experimental  current  for  the  second  test. 
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Float  Trapping  Efficiency  Modal.  This  nodal  consider*  the  flux 
trapping  efficiency  without  regard  to  the  flux  loas  mechanism.  A  plot 
of  the  efficiency  for  the  first  test  ia  shown  in  Figure  6-10  and 
similar  curves  for  the  remaining  testa  ere  in  Appendix  B.  Points  1» 

2,  end  3  from  Figure  3-4,  corresponding  to  key  points  in  generator 
operation,  are  shown  in  this  figure. 

In  each  teat,  the  efficiency  remains  better  then  90S  for  the  firat 
half  of  generator  operation.  The  values  greater  than  100Z  in  the  firat 
experiment  again  indicate  that  flux  compression  ia  occurring  faster 
than  that  predicted  by  the  Gurney  approximations.  During  the  next  25X 
of  generator  operation  the  efficiency  drops  to  approximately  80X. 

e 

During  the  final  20  microseconds,  severe  losses  occur.  The  final 
efficiency  ranges  between  50  and  60  percent. 

The  curves  from  the  first  three  experiments  were  used  to  generate 
a  best  and  worst  case  efficiency  model.  The  plots  for  the  best  and 
worst  case  currents  for  the  efficiency  model  for  the  first  test  are 
shown  in  Figure  6-11.  Curves  for  the  remaining  experiments  are 
presented  in  Appendix  B.  At  the  end  of  generetion,  the  final  currents 
in  each  test  fall  well  within  the  predicted  range  indicating  that  the 
efficiency  model  ia  reliable. 

General  Results.  All  the  models  reliably  predict  the  generator 
output  for  the  firat  50Z  of  generator  operation.  Beyond  that  point 
only  the  empirical  models  can  predict  the  current  values.  Baaed  on  the 
experimental  data,  the  computer  model  can  reliably  predict  the  range 
of  output  currents. 
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Figure  6-11.  Efficiency  Currents  Test  1 
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Conclusions 

Several  conclusions  can  be  drawn: 

1)  The  computer  model  can  predict  the  generator  output  current 
reliably  during  the  first  50%  of  operation,  it  can  also  predict  the 
ranges  in  which  the  final  maximum  output  current  will  fall* 

2)  The  losses  during  the  initial  SOX  cf  generator  operation  can 
be  neglected  and  the  generator  may  be  considered  lossless  for  this 
period  J 

3)  The  plate  resistance  and  flux  diffusion  resistance  account 
for  only  a  minor  portion  of  the  loss  mechanisms  in  generator  operation; 

4)  The  empirical  models  show  that  the  flux  compression  described 
by  the  Gurney  approximations  is  net  totally  accurate; 

5)  Losses  are  severe  during  the  final  moments  of  generator 
operation;  and 

6)  The  validation  experiments  show  that  the  strip  MCG  is  capable 
of  delivering  large  currents  over  a  relatively  long  period  of  time. 

Recyrnggindat  ions 

Rased  on  the  work  performed  for  this  thesis  several  recommendations 
can  be  made.  These  include: 

1)  Computet  validation  of  the  trapezoidal  portion  of  the  computer 
program  has  never  been  accomplished.  This  portion  of  the  program 
should  be  validated; 

2)  The  empirical  models  are  based  on  the  results  of  only  three 
teats.  Further  tests  should  be  conducted  to  further  validate  them  or 
to  refine  their  accuracy; 
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3)  Strip  generators  other  than  the  AFWL  MCG  should  be  used  to 
further  validate  the  models  (for  example ,  different  materials  or 
dimensions) ; 

4)  Because  the  loss  mechanisms  in  bellows  generators  are  similar , 
the  computer  program  should  be  expanded  to  include  this  geometry  MCG; 

5)  Because  the  generator  drives  a  low  inductance  load,  the  use 
of  transformers  to  drive  capacitive  and  resistive  loads  should  be 
investigated  and  perhaps  incorporated  into  the  computer  program; 

6)  Due  to  the  momentum  balances  performed,  tapered  explosives 
should  be  considered  to  overcome  the  intense  magnetic  pressure  as 
generator  operation  draws  to  an  end; 

7)  When  high  magnetic  fields  and  high  currents  were  present,  the 
load  coil  was  visibly  deformed;  it  opened  up.  This  effect  should  be 
included  in  the  computer  program; 

8)  High  speed  photographs  of  generator  operation  should  be  done 
to  gain  a  better  understanding  of  how  the  inductance  actually  changes 
with  time; 

9)  Explosive  type  and  plate  material  should  also  be  varied  to 
further  validate  the  computer  model;  and 

10)  Due  to  the  intense  magnetic  pressures  present  near  the  end  of 
generator  operation,  a  momentum  balance  should  be  incorporated  into 
the  computer  program. 


Appendix  A 


This  appendix  describes  the  computer  program*  This  program  is 
vritteu  in  FORTRAN  V.  The  input*  to  the  program  are  best  described 
by  Figure  A-l. 


Inputs 

TPFW  is  the  parallel  widch  of  the  top  plate.  TPTW  is  the 
trapezoidal  width  of  the  top  plate  at  the  edge.  TPT  is  the  plate 
thickness.  XPTL  is  the  trapesoidal  length  of  the  top  plate.  TPL  is 
the  total  length  of  the  top  plate.  TTH  is  the  height  of  the  trough 
of  the  top  plate.  BPL  is  the  length  of  the  bottom  plate.  CWH  is  the 
height  of  the  cutting  wedge.  CWTL  is  the  length  of  the  top  of  the 
cutting  wedge.  CWBL  is  the  length  of  the  bottom  of  the  cutting  wedge. 
JWTL  is  .-he  length  of  the  top  of  the  jet  wedge.  JWBL  is  the  length  of 
the  bottom  of  the  jet  wedge.  JCL  is  the  length  of  the  jet  catcher. 

THT  is  the  throat  thickness  at  the  jet  catcher.  TTH  is  the  width 
of  the  steel  tamper.  TTT  is  the  thickness  of  the  top  steel  tamper. 

TTL  is  the  length  of  the  steel  taaiper.  EXPL  is  the  weight  of  the 
explosive.  CEMW  is  the  weight  of  the  cement  brick  tampers.  VOLT  is 
the  capacitor  bank  voltage.  LINEL  ia  the  line  inductance.  LOADL  is 
the  load  inductance. 

The  load  geometries  are  required  for  the  path  length  calculations 
They  can  be  found  using  Figure  A-2  as  a  guide.  RLOAD  is  the  radius  of 
the  load  coil.  For  small  loads  RLOAD  is  zero.  DLOAD  is  the  depth  of 
the  load.  WLOAD  is  the  width  of  the  load. 

Two  filenames  are  also  inputs.  They  are  the  best  case  model 
results ,  FILNAM  and  the  worst  case  model  results,  FNAME. 


Outputs 

Outputs  for  this  program  consist  primarily  of  currents  and 
fields.  TAMH  is  the  total  tamper  mass.  NTOC  is  the  metal  mass  to 


explosive  ratio.  NTOC  is  the  tamper  mass  to  explosive  ratio.  The 
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Figure  A-2 .  Load  Geometries 


detonation  velocity  is  D.  The  metal  velocity  is  VA.  The  Gurney  angle 
is  THETD.  TIMEL(x)  is  the  time  varying  Gurney  inductance.  TIMEI(x) 
is  the  flux  conserved  current.  HAGH(x)  and  MAGB(x)  are  the  magnetic 
field  and  the  magnetic  flux  density  respectively  for  the  flux  conserved 
model.  RES(x)  is  the  time  varying  resistance.  IRES(x)  is  the  current 
for  the  resistance  model.  KAGRH(x)  and  MAGRB(x)  are  the  magnetic 
field  and  magnetic  flux  density  respectively  for  the  resistance  model. 
ISLIP(x)  is  the  current  for  the  sliding  contact  loss  model.  HSLIP(x) 
and  BSLIP(x)  are  the  magnetic  field  and  the  magnetic  flux  density 
respectively  for  the  sliding  contact  loss  model.  lEFF(x)  is  the 
efficiency  model  current.  HEFF(x)  and  BEFF(x)  are  the  magnetic  field 


and  magnetic  flux  density  respectively  for  the  efficiency  model. 

IAA(x)  is  the  changing  area  model  current.  SAA(x)  and  HAA(x)  are  the 
magnetic  flux  density  and  magnetic  field  respectively  for  the  changing 
area  model. 

The  program  prints  tabular  output  files  of  these  variables.  It 
additionally  creates  eight  output  files  specifically  for  plotting  the 
currents  described  in  Chapter  VI.  A  complete  copy  of  the  computer 
program  can  be  found  on  the  following  pages. 
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* 

* 

# 

*  this  program  calculates  and  prints  the  ejected  results  of  a  strip 

*  GENERATOR  GIVEN  THE  BANK  CAPACITANCE.  LINE  INDUCTANCE,  LOAD  INDUCTANCE. 

*  GENERATOR  DIMENSIONS.  EWLOSIVE  SIZE.  AND  TAMPER  SIZE.  AT  PRESENT, 

*  THIS  PROGRAM  ASSUPES  THAT  THE  GENERATOR  IS  MADE  OF  COPPER  AND  THE 

*  EXPLOSIVE  IS  DETASHEET.  IF  YOU  UISH  TO  CHANGE  TO  ANOTHER  MATERIAL  OR 

*  EXPLOSIVE,  PARAMETERS  WHICH  MUST  BE  CHANGED  ARE  INDICATED  BY  COMMENTS 

*  SURROUNDED  BY  A  ROW  OF  *'S.  THIS  PROGRAM  UTILIZES  THE  GURNEY  EQUA- 

*  TICNS  AND  SIPPLE  GEOMETRY.  ALTHOUGH  THE  PROGRAM  ALLQUF  FOR  A  TRHPE- 

*  ZOIDAL  GENERATOR,  IT  ASSUMES  THAT  THE  STEEL  TAPPER  IS  RECTANGULAR  IN 

*  CROSS  SECTION.  THIS  PROGRAM  CALCULATES  BOTH  A  FLUX  CONSERVED  MODEL 

*  AHD  A  TIME  VARYING  RESISTANCE/INDUCTANCE  MODEL  BASED  ON  THE  MAGNETIC 

*  SKIN  DEPTH.  THIS  PROGRAM  ALSO  CALCULATES  THREE  EPPIRICAL  MODELS.  THE 

*  SLIPPING  LOSS  EPPIRICAL  MODEL,  THE  FLUX  TRAPPING  EFFICIENCY  MODEL,  AND 

*  THE  CHANGING  AREA  EPPIRICAL  MODEL. 

* 

* 

PROGRAM  STRIP 

INTEGER  LIMIT.  L1M1.LIMI2.LIM2.THLIH,  I, J.K.L.M.N,P.Q,R,S.T.U.V 
INTEGER  II,  JJ.KK.LL.Pli.  IJ,  IK.CK1.CK2 

REAL  TPU.TPT.TPL.TTH.BPL.CUH.CUTL.CIJBL, JUTL, JUBL,JCL.TnJ,TTT,TTL 
REAL  THT.CEMU. EXPL. VOLT. CAP. L INEL. LOADL. VOL  1 . VOLZ. PVOL,  PETU. TAPMJ 
REAL  TAMU.  HTQC .  NTOC . AHUM. ADEN, A. VMI , VM2.  VM.  VNUM.  VDEN.  THETA.  ARC  t ,  D 
REAL  XR.Xl.XE.XFt,  INTL.THROL.XSTOP.TSTOP.TPTL.TPTU.RESTl.CUBM.AA 
SEAL  ATQT.  ARG2.  APR  IM.  THETD.  ALPHD.  PH  ID.  PH  I ,  ZETA,  ZETD .  ZETPD.  ZETAP  ,  X 
REAL  LALPH. ALPHA. Y1 . ASUB I . AGEN. XQ. GAMAD. GAMA. ASUB2. ASUB3 .LGAPIA, VA 
REAL  XDELT, H3. CENT1 . CENT2, CENT , CHUM, CDEN,  IQ.TOTL.  IT.TPPU.BASE 
REAL  SBASE, SHT. TRAPA.  IMHT.ZREF.CALH.CALU.>«_IM2.XLlrB.Bl,B2,Hl.H2 
REAL  INTL1.  INTL2,  INTL3.  INTL4.ANGLE.AI.LHT,V0L4.L1,L2.XLIN4,SLA1 
REAL  SLA2. SLA3. SLA4.SLA5.SLA6.  H.IM5.H.  IMl . TPPLM, LA  1 . RLOAD, DLOAD 
REAL  ULOAD .  D IF . LPH I .  LZET.  TLOAD.  PMAG,  ASUB4.  ASUB5,  MU .  REST,  L I TK, BB 
REAL  REST2.SL I FD, POWER. CONS. ROEB. CC. DD. EE, FF. GG, DEC  I . DEC2,  DEC3 
REAL  DEC4.DEC5, DECS, DEC7.D1.D2.D3.D4. LEX 

REAL  AREA (SB8 )  .MISEC  (586) .  TIMEL  (568) .  TIME  I  (500) ,  UC500) ,  LA4(508) 
REAL  LAS (580) , LAS ( 500) . LA4D ( 500) . LASD ( 508) , LA6D ( 580 ) .  POS  (500) 

REAL  LA7  ( 500) .  LA7D  ( 500) .  ARE  A 1  <  508)  .  AREA2  (500) .  AREA3  ( 500 ) .  TAU  (500) 
REAL  ARE  ID  (508)  ,  ARE2D  (500) ,  ARE3D  (588) ,  TVL 1  (500) ,  TVL2  (580) ,  KO  (500) 
REAL  AREA4(508) . AREA5(500) , ARE4D(5B6) . ARE5DG08) ,TVL4(50B) 

REAL  TVL3(500),TLeN(50B>,RLEN(500),LLEN(508).BLEN(50B).SlGMA(500) 
REAL  PATH (500) .  MAGB (500) ,  PIAGH (500) , TVL5 (500) , RES (500) ,  TEMP (508) 
REAL  RAV(  500) ,  L AV  ( 500 ) ,  EARG  1(508),  EARG2  ( 500 ) ,  I  RES  ( 500 ) ,  MAGRH  ( 500 ) 
REAL  MAGRB  (500) .  EARG3  (500) ,  LA8  (580)  ,LA8D(500) ,  EXPO  (500) .  ROE  (500) 
REAL  SMS (500) .EFF (500) . AAA (580) , ISL IP (500) , IEFF (500) ,  I AA (580) 

REAL  HSL  IP  (500)  .  BSL  IP  (500) ,  HEFF  (500) ,  BEFF  (588) .  HAA  (580) ,  BAA  (500) 
REAL  SPHK500) 

CHARACTERS  FNAPE 
CHARACTERS  F  ILHAM 
CHARACTERS  PLOT1 
CHARACTERS  PLQT2 
CHARACTERS  PLOT3 
CHARACTERS  PLOT4 
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CHARACTER*^  PLOTS 

* 

*  THIS  PORTION  OF  THE  PROGRAM  HAS  THE  USER  ENTER  THE  NEEDED  PARAMETERS. 


PRINT*, 'PLEASE  ENTER  THE  WIDTH  OF  THE  PARALLEL  PORTION  OF  THE' 
PRINT*, 'TOP  PLATE  IN  CENTIMETERS. >' 

REAT*,TPPU 

PRINT*, 'PLEASE  ENTER  THE  WIDTH  OF  THE  TOP  PLATE  AT  THE' 

PRINT*, 'TRAPEZOIDAL  END  IN  CENTIMETERS. >' 

READ*,  TPTU 

PRINT*, 'PLEASE  ENTER  THE  THICKNESS  OF  THE  TOP  PLATE' 

PRINT*,' IN  CENTIMETERS. >' 

READ*,  TPT 

PRINT*, 'PLEASE  ENTER  THE  LENGTH  OF  THE  TRAPEZOIDAL  PORTION  OF  THE' 
PRINT*, 'TOP  PLATE  IN  CENTIMETERS. >' 

READ*, TPTL 

PR iNT*, 'PLEASE  ENTER  THE  TOTAL  LENGTH  OF  THE  TOP  PLATE' 

PR INT*, ' IN  CENTIMETERS . >' 

READ*  TPL 

PR INT*, 'PLEASE  ENTER  THE  HEIGHT  OF  THE  TROUGH  IN  CENTIMETERS. >' 
READ*,TTH 

PR INT*, 'PLEASE  ENTER  THE  LENGTH  OF  THE  BOTTOM  PLATE' 

PRINT*,' IN  CENTIMETERS. >' 

READ*,BPL 

PRINT*, 'PLEASE  ENTER  THE  HEIGHT  OF  THE  CUTTING  WEDGE' 

PRINT*,' IN  CENTIMETERS. >' 

READ*,CWH 

PR INT*, 'PLEASE  ENTER  THE  LENGTH  OF  THE  TOP  BASE  OF  THE  CUTTING' 
PRINT*, 'WEDGE  IN  CENTIMETERS. >' 

READ*,CUTL 

PR INT*. 'PLEASE  ENTER  THE  LENGTH  OF  THE  BOTTOM  BASE  OF  THE  CUTTING' 
PRINT*. 'WEDGE  IN  CENTIMETERS. >' 

READ*  CUB, 

PR INT*. 'PLEASE  ENTER  THE  LENGTH  OF  THE  TOP  BASE  OF  THE  JET  WEDGE' 
PRINT*,' IN  CENTIMETERS. >' 

READ*, JUTL 

PR  I NT*. 'PLEASE  ENTER  THE  LENGTH  OF  THE  BOTTOM  BASE  OF  THE  JET' 
PRINT*. 'WEDGE  IN  CENTIMETERS. >' 

READ*,JWBL 

PR INT*. 'PLEASE  ENTER  THE  LENGTH  OF  THE  JET  CATCHER' 

PRINT*,' IN  CENTIMETERS. >' 

READ*. JCL 

PRINT*, 'PLEASE  ENTER  THE  THROAT  THICKNESS  AT  THE  JET  CATCHER  IN' 
PRINT*, ' CENTIMETERS . >' 

READ*,THT 

PRINT*, 'PLEASE  ENTER  THE  WIDTH  OF  THE  TOP  STEEL  TAMPER' 

PRINT*, ' IN  CENTIMETERS . >' 

READ*,TTW 

PRINT*. 'PLEASE  ENTER  THE  THICKNESS  OF  THE  TOP  STEEL  TAMPER  IN' 
PRINT*, 'CENTIMETERS. >' 

READ*, TTT 

PR INT*. 'PLEASE  ENTER  THE  LENGTH  OF  THE  TOP  STEEL  TAMPER' 

PRINT*. 'IN  CENTIMETERS. >' 

READ*.TTL 


PRINT#, 'PLEASE  ENTER  THE  WEIGHT  OF  THE  E>PLGSIVE  IN  POUNDS. >' 
READ*,EXPL 

PRINT#. 'PLEASE  ENTER  THE  WEIGHT  OF  THE  TOP  CEMENT  BRICKS  IN* 
PRINT*. 'POUNDS.)* 

READ#,CEMU 

PRINT*. 'PLEASE  ENTER  THE  BANK  VOLTAGE  IN  VOLTS. >' 

READ*, VOLT 

PRINT*. 'PLEASE  ENTER  THE  BANK  CAPACITANCE  IN  MICROFARADS. >' 

READ*. CAP 
CAP -CAP* IE- 06 

PRINT*. 'PLEASE  ENTER  THE  LINE  INDUCTANCE  IN  NANOHENRIES. >' 
READ*,L!NEL 

PRINT*. 'PLEASE  ENTER  THE  LOAD  INDUCTANCE  IN  NANOHENRIES.)' 
READ*,LOADL 

PRINT*. 'PLEASE  ENTER  THE  RADIUS  OF  THE  LOAD  IN  CENTIMETERS.)* 
READ*.RLOAD 

PRINT*,' PLEASE  ENTER  THE  DEPTH  OF  THE  LOAD  IN  CENTIMETERS.)' 
READ*,DLOAD 

PRINT*. 'PLEASE  ENTER  THE  WIDTH  Or  THE  LOAD  IN  CENTIMETERS.)* 

READ*,WLOAD 

LINEL-LINEL*IE-09 

LOADL ■LOADL* IE-09 

PRINT*. 'PLEASE  ENTER  THE  BEST  RESULTS  FILENAME.)* 

READ  (5.500)  FILNAM 
500  FORMAT  (A48) 

PRINT*. 'PLEASE  ENTER  THE  WORST  RESULTS  FILENAME.)* 

READ  (5.500)  FNAME 

* 

*  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  MASS  AND  TAMPER  RATIOS  AND 

*  USES  THE  GURNEY  EQUATIONS  TO  DETERMINE  A  METAL  VELOCITY. 

* 

VOL * -TPPU*(TPL-TPTL) *TPT 
BASE  "TPPU/2 . 0 
SBASE-TPTU/2.0 
SHT-BASE-SBASE 
IF  (SHT.EQ.0.0)  THEN 
VOL2-0.0 
ELSE 

TRAPA-ATAN (SHT/TPTL) 

I MHT-8  ASE^TAN  ( TRAPA ) 

LHT-TPTL-IMHT 

AI«IMHT*8ASE 

A2-LHT*SBASE 

VOL2-TPT*(Al-A2) 

END  IF 

Ll-SQRT(  ( IMHT*IMHT)  +(BASE>«ASE) )  -SQRT(  (LHT*LHT)  +(SBASE*SBASE) ) 
L2-TPL-TPTL 

V0L3-2 . 0*TTH*(L  1+L2)  *TPT 
P VOL "VOL 1+VOL2+VOL3 
METU-  (B .  96*P  VOL)/' 1000 . 0 
CEMU-CEMU/2.  20462 
V0L4-TTU*TTT*TTL 
TAMMJ-  (7 . 86*V0L4)  /1000 . 0 
TAMW-TAMIMILHCEMU 
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‘  C  V*  V 


A  •* 


EXPL-EXPL'2. 20462 
MTOC-r-ETU'EXPL 
NTOC  -TAMU/EXPL 
ANUM-1.0+(2.0*MTOC) 

ADEN-1.0+(2.0*NTOC) 

A-ANUM/ADEN 

vm-u.0+(fl*«#«))/(3.0*ci.0+fl)) 

VM2-NT0C#A#A 

VM3-MTOC 

*  VNUM  WILL  HAVE  TO  BE  CHANGED  IF  THE  Etf>LOSIVE  IS  NOT  DETASHEET. 


VNUM-2.2? 

VDEN-SQRT(VM1+VM2+VM3) 

VM-VNUIWDEN 

* 

*  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THETA  AND  V  USING  KENNEDY'S 

*  fCTHOD  AND  THE  GURNEY  EQUATIONS. 

* 


ilaMBMololoiQicWnto^ 

*  D  UILL  HAVE  TO  BE  CHANGED  IF  THE  EXPLOSIVE  IS  NOT  DETASHEET. 
iMnjnMotolnWnMnioMranlnMnMninMoMnKiMnlolnloMBWoM^ 


>  D-6.B 

ARG1-VTV(2.0*D) 

THETA-2 . 0*<  AS IN (ARG l ) ) 

VA-D*TAN (THETA) 

* 

i  *  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  CHANGING  WIDTH  OF  THE 
j  *  GENERATOR. 

* 

XE-BPL-CUTL-JCL 

XE-XE*1E81 

XSTOP-XE+C  (CUH*1E0  H/TANOHETA) ) 

TSTOP-XSTOP/D 
LIMIT-TSTOP+1 
,  LIM2-2*LIMIT 

I  TUt-0.5 

DO  5  l-l, LIMIT, I 
X-D*TIME 
POS(I)-X 
XL  IM1  -TPTL*lE0  1 
IF  (X.LE.XLIM1)  THEN 
BASE-TPPU/'2.0 
SBASE-TPTU/2 . 0 
SHT-BASE-SBASE 
TRAPA-ATAN (SHT/TPTU 
I MMT -BASE/TAN ( TR APA ) 

LHT  -  (TPTL- IMHT) #1E0 1 
CALU- CLHT+X) *TAN ( TRAP A) 

U(I)-CALU*lE-03*2.0 

ELSE 

U(I)»TPPW*lE-02 
END  IF 

TIfE-TirE+1.0 
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S  CONTINUE 

m 

*  THIS  portion  of  the  program  calculates  the  initial  inductance  and  the 

*  THROAT  inductance. 

* 

XR-CUBL-CUTL 
TPTLH-TPTL*1E01 
CUBM-CUBL*1E01 
ALPHA -ATANtXR/ClJH) 

N_IM2-(TPL-JCL)*1E01 
XLIM3-(TPL-JU9L)*1E01 
XLIM4»<TPL-JUTL)*1E01 
DO  10  J-l. LIMIT, l 

IF  (TPTL.EQ.0.0)  THEN 
SLA l -0 . 5*XR*CUH* IE-04 
Xl-BPL-CUBL-JCL 
SLA2-CLJH*Xl*lE-04 
SLA3-CUH*(JCL-JLJBL)#1E”04 
XF1-JUBL-JUTL 

SLA4-0 . 5*XFl*(CUH-TPT)*tE-04 
SLA5-SLA 1+SLA2+SLA3+SLA4 
SLA6-JCL*THT*lE-04 
ELSE 

XR-XR*tE01 
IF  (X.LE.XR)  THEN 

LA4(K)-0.5*X*(X/TAN(ALPHA).>*1E-0G 
LA5CK) -0.0 
LAS  (0-0.0 
LA7 00-0.0 

ELSE  IF  (X.LE.XLIM2)  THEN 

LA4(K)-0.5*XR*<>R/TAN<ALPHA))*lE-0b 
LAS  OO  • fX-XR)  *CUHM*iE-06 
LA6 00-0.0 
LA70O -0.0 

ELSE  IF  (X.LE.XLIM3)  THEN 

LA4 00  >0 . 5*XR* ( XR/TAN ( ALPHA ) ) * I E-06 
LA50O -(XL  IM2-XR)  *OJHM*lE-06 
LA6  OO  -  < X-XL IM2) *(  (CUH-THT-TPT) *1E01 ) *IE-06 
LA7(K) -0.0 
ELSE 

ANGLE -ATAN <  UIBL-JUTL)/CUH) 

H1-CX-XLIM3)*1E-01 
B1 -CUH-THT-TPT 
H2-JUBL-JUTL-H1 
B2  -H2/'  ( TAN  ( ANGLE ) ) 

LA40O  -0 . 5*H?*(>«/TAN  (ALPHA) )  *lE-06 
L.A50O -(XL  IM2-XR)  *CUHM*tE-06 
LA6  (K )  -  04.  IM3-H.  IM2 )  *<  (CUH-THT-TFT)  *1E0 1 )  *1E-0S 
LA7OO-0.5*CBl+B2)*Hl*lE-04 
END  IF 
XR-XR*lE-01 
END  IF 

TlME-TIHE+1 .0 
10  CONTINUE 
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DO  15  K-2.LIMIT.1 
L-PHl 

LA4D(1)-LA4(1) 

LA5DU)-LA5<1) 

LA6D(1)-LA6(1) 

LA7D<1)-LA7U) 

LA4DCK) -LA4<K)-LA4(L) 
LA5D(K)-LA5<K>-LA5(L) 
LA6DCK)  -LA6(K)-LA6(L) 


LA7D<K)«LA7<K)-LA7<L> 

15  CONTINUE 

DO  20  M-t.LIMIT.l 

INTL  l  -  ( ( l .  256637 1*LA4D  CM)  >  /UtM) )  MNTL 1 
I NTL2  -  ( <  1 . 256637 1  *L  A5D  (M) )  /UtM) ) + 1 NTL2 
I NTL3- ( C  l .  256637 1*L  A6D  (II) )  /UCM) )  +  INTL3 
INTL4-((  1 .256637 1*.A7D  CM) ) /UCM)  )+INTL4 
20  CONTINUE 

IF  (TPTL.EQ.0.0)  THEN 

INTL-(1.2566371*SLA5)/CTPPU*lE-02) 

INTL- INTL* IE -06 
ELSE 

INTL- INTL 1+INTL2+INTL3+INTL4 
INTL- INTL*lE-06 
END  IF 

THL IM-( CTPL-JCL) *1E0 1 )  /D 
TIME -0.0 

DO  25  N-l. LIMIT. 1 
XLIM5-THLIP*D 
X-D*TII€ 

IF  CX.LE.XLIM5)  THEN 
LA0(N) -0.0 
ELSE 

LA8 (N) - ( ( CBPL-JUBL)*1E0 1 ) -X) *CTHT*1E0 1 ) *lE-06 
END  IF 

TIIt-TIME+1.0 
25  CONTINUE 

DO  30  P-2. LIMIT, 1 
Q-P-l 

LA8DC  D-LA8C1) 

LA8DCP) -LAS (P) -LAB (Q) 

30  CONTINUE 

IF  (TPTL.EQ.8.0)  THEN 

THROL-( l . 256637 l*SLA6)/(TPPW*lE-02) 

ELSE 

DO  35  R-l. LIMIT, 1 

THROL  ■  (1 . 256637 1  *.  A8D  (R) )  AH  R)  +THROL 
35  CONTINUE 

END  IF 

THROL -THROL* IE-06 

*  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  HOU  THE  INDUCTANCE  CHANGES 

*  UITH  TIME.  IT  ASSUrES  A  CONSTANT  DETONATION  VELOCITY.  IT  FURTHER 

*  ASSUPES  THAT  THE  PETAL  LINER  IS  TOTALLY  ELASTIC  AND  THAT  THE  LOAD  AND 

*  THROAT  INDUCTANCES  DO  NOT  VARY  UITH  TIPE. 
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XE-BPL-CUTL-JCL 

XE-XE*IE01 

XSTOP»XE+(C(ClJH-THT-TPT)*lE0t) /TAN  (THETA)) 

TSTOP-XS TOP/D 

ARG2«CUH/(CUBL-CUTL> 

ALPHA -ATAN(ARG2) 

THETD«THETA*(l80.0/3. 141592?) 

ALPHD-ALPHA*( 180.0/3. 1415927) 
PHID*180.0-THETD-ALPHD 
?HI-PHID*(Z. 1415928/180.0) 

ZETrt-ATAN<XFl/CUH) 

ZETD»Z£TA*( 180.0/3. 1415927) 

2ETPD-90 . 0+ZETD 
ZETAP-ZETPD*(3 . 1415927/180.0) 

ALPHA-ALPHD*<3. 1415927/180.0) 

XR-CUBL-CUTl 
XR»XR*1E01 
Xl-BPL-CUBL-JCL 
XF 1  ■ JWBL-JUTL 
XO-CUH/TAN (THETA) 

XO»XO*1E01 

LimT-Tsrop+i 

TIME -0.0 

DO  40  S-l, LIMIT, 1 
X»D*TIME 
XLIM-XR+XO 

XF  -  (BPL-CUTL- JUTL ) * 1 E0 1 
IF  (TPTL.EQ.S.0)  THEN 

IF  Os.LE.XLIM)  THEN 
LALPH-SIN(THETA)*(X/(SIN(PHI) ) ) 

Y 1 -LALPH*S I N  <  ALPHA) 

ASU81-0.5*X*Y1 
AGEN- <SLA5*1E06) -ASUB 1 
AREA(S) -AGEN*lE-06 
MISEC(S) -TIME 
ELSE  IF  (X.LE.XF)  THEN 
ASUB 1-0.5**.  IM*CUH*1E0 1 
ASUB2 » ( X-XL I M) *CUH* ) E0 1 
AGEN -<SLP5*1E06> -ASUB l -ASUB2 
AREA(S) «AGEN*1E-0S 
MlSEC(S) -TIME 
ELSE 

ASUB  1  -0 . 5**.  I  M*CUH*1E0 1 
ASUB2-CX-XLIM) *CUH*1E01 
GAMAD -180. 0-ZETPD-THETD 
GAMA-GAMAD*(3. 1415927/180.0) 
XDELT-X-XF 

LGAMA-S IN (THETA) *(XDELT/S IN(GAMA) ) 
H3-LGPMA*S IN (GAMA) 

:NTI -SQRT< (LGAMA*LGAMA) - (H3*H3) ) 
CENT2-SQRT( (XDELT*XDELT)-(H3*H3) ) 
CENT -CENT 1 +CENT2 
ASUB3 -0 . 5*CENT*H3 


TSUB4-  '  (THT*TPT)  *<  JCI.-J:J3L)  *lE02) 

A3U85  » . 5*(TPT*TPT* 1 E32 ) /TAN  < THETA) 
AGEN»(5LA5*1E0S)-ASUB 1-ASUB2+ASUB3-ASUB4+ASUB5 
AREA(S)  =AGEN*tE-06 
MI3ECCS) “TIME 
END  IF 

TIME-TIME+1.0 

ELSE 

IF  (X.LE.XLIM)  THEN 

LALPH-SIN(THETA)*(X/(STN(PHI))) 

Y1-LALPH*S IN (ALPHA) 

ASU81«0.S*X*Y1 
AREA  1 <S) *ASUB 1*1E”06 
AREA2(S) *0.0 
AREA3(S) *0.0 
AREA4(S) *0.0 
AREAS (S) >0.0 
MISEC(S) »TIME 
ELSE  IF  (X.LE.XF)  THEN 

AREA  1 (S) -0 . 5*XL IM#CUH*1E0 1*1E-0G 
AREA2 (S) - <  X-XL I M) *CUH* 1E0 1 *lE-06 
AREA3(S) "0.0 
AREA4CS) *0.0 
AREAS (S) >0.0 
MISEC(S) *TIME 
ELSE 

AREA1  (S)  >0 . 5*XL  IM*CUH*1E0  l*lE-06 
AREA2CS)  ■<X“XLH1)*CUH#1E01*1E“06 
GAMAD- 180 . 0-ZE7PD-THETD 
GAMA-GAf1AD*(3. 1415927/100.0) 

XDELT-X-XF 

LGAMA-S !N (THETA) *(XDELT/S IN (GAMA) ) 
H3*LGAMA*SIN(GAMA) 

CENT1 «SQRT( (LGAMA*LGAMA) - (H3#H3> ) 

CENT2-SQRK  ( XDELT*XDELT)  -  ( H3*H3 ) ) 
CENT-CENT1+CENT2 
AREA3 (S) >0 . 5*CENT*H3*1E-0S 
AREA4 ( S ) - ( ( THT*TPT) *(JCL-JUBL))*1 E-04 
AREAS (S) - (0 . 5*(TPT*TPT> /TAN (THETA) ) #lE-04 
MISEC(S)»TIME 
END  IF 

TIME“TIME+1 .0 
END  IF 
40  CONTINUE 

IF  (TPTL.EQ.0.0)  THEN 
DO  45  T-l, LIMIT. 1 

TVL*((l. 256637 1*AREA <T) ) /(TPPU*lE-82) ) *lE-06 
TIMEL(T) -TVL 
45  CONTINUE 

ELSE 

DO  50  U"2. LIMIT,  1 
V-U-l 

ARE  ID ( 1) -AREA1 (1) 

ARE2D( 1) »AREA2( 1) 
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ARE3D(1)-AREA3(1) 

ARE4D  < 1 ) -AREA4( 1 ) 

ARE5DC  t) -AREA5( 1) 

ARE  ID (U) -AREA  1 (U) -AREA  1 (V) 

ARE2D (U) -AREA2 <U) -AREA2 C V) 

ARE3DCU) -AREA3  <U) -AREA3 ( V) 

ARE4D  OJ ) • ARE A4 < U ) - ARE A4 ( V) 

ARE5D (U) -AREA5 (U) -AREAS (V) 

50  CONTINUE 

DO  55  1-2, LIMIT.  1 
J-I-l 

TVL1 ( 1)  ■(  ( 1.256637 1#ARE1D( 1) )/U( 1) )*1E-B6 
TVL2 ( 1) • ( ( 1 . 256637 l*ARE2D (1 ) )  Al( 1 ) ) *lE-06 
TVL3U)-<(1. 256637 1*ARE3DU) ) A1C1)  )*lE-06 
TVL4< 1 ) - < (1 . 256637 1*ARE4D  U ) )  Al( 1) ) *lE-06 
TVL5 (1 ) - < (1 . 256637 1*ARE5D (1) )  AH  1 ) ) *IE-06 
TVL 1 ( I) - ( ( (1 . 256637 1 *ARE ID ( I ) ) /U( I) ) * 1 E-06) +TVL 1 ( J) 

T VL2 <  I ) ■ < < < 1 . 256637 1 *ARE2D (I))/U(I))*1 E-06) +TVL2 < J ) 

TVL3 ( I ) - ( ( < 1 . 256637 1  *ARE3D C I ) ) /UC I ) ) ♦ 1 E-06) +TVL3 ( J ) 

,  TVL4<  I ) « <  <  < 1 . 256637 1 *ARE4D ( I ) )  /U( I ) ) * I E-06 ) +TVL4( J) 

TVL5 ( I) - C  < ( 1 . 256637 1*ARE5D ( I ) ) /U( I ) ) * 1 E-06) +TVL5 ( J) 

55  CONTINUE 

DO  60  K-l. LIMIT,  l 

TVL- INTL-TVL 1 (K) -TVL2CK) +TVL3 (K) +TVL4(K) -TVL5(K) 

TIMELCK) -TVL 
60  CONTINUE 

END  IF 

* 

*  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  AND  PRINTS  THE  CHANGING  CURRENT 

* 

DO  65  L-l, LIMIT, 1 

CNUM- INTL+THROL+LOADL 
CDEN-THROL+LQADL+TIMEL(L) 

TOTL-CNLM+L INEL 
IO-VOLT*(SQRT(CAP/TOTL) ) 

I T- ( CNUM/CDEN) * I 0* 1 E-03 
TIMEI(L) -IT 
TIMEL(L) -TIMEL(L)*IE09 
65  CONTINUE 
♦ 

♦THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  CHANGING  PATH  LENGTH  OF 

♦  THE  SYSTEM. 

♦ 

IF  (RLOAD.EQ.0.0)  THEN 

TLOAD-  (2 . 0^3 . 1415327*<RL0AD) +( (2 .  O^DLOAD)  HLOAD)  1 
ELSE 

TLOAD- (2 . 0*<3 . 1415927*<RL0AD)  -UL0AD+I2 .  O^DLOAD) 

END  IF 
TIME-0.0 

DO  70  L-l. LIMIT, 1 
X-D^TIfE 

IF  (X.LE.N.IM)  THEN 

TLEN  CL)  -TPL-JUTL-CUTL-  CX*lE-0 1 ) 
LALPH-SIN<THETA)^(X/(SIN(PHI)miE-0l 


DIF-CLBL-C'JTL 

LPH  X  aS IN l ALPHA)  *(X/(S IN (PH  I ) ) ) *lE-9 1 
LLEN(L)  -(SORTI  (CUH*CLJH)+(D  IF*DIF) )  )-LALPH+LPHI 
BLEN(L) -TPL-CUBL+TLOAD 

RLEN(L)  -SGRT(  ( ( JUBL-JUTL)  *(  JUBL-JUTL) )  +(CUH*CLJH) )  +JCL 
ELSE  IF  (X.LE.XF)  THEN 

TLEN(L) -TPL-JUTL-CUTL- <X*lE-0 1 ) 

LLEN(L) -CUH/S IN (THETA) 

BLEN  i.L)  -TPL -CUTL+(CUH/TAN ( THETA) )  +TL0AD- CX*1E-01) 

RLEN(L)  -SQRT(  ( <  JLJBL-JUTL)*(  JUBL-JUTL) )  +(CUH*CUH) )+JCL 
ELSE 

TLEN(L) “0.0 
XDELT»(X-XF)*IE-01 
GAMAD- 180 .0-2ETPD-THETD 
GAMA-GAMAD*<3.  141592?/'180.0) 

L2ET«SIN(2ETAP) *(XDEIT/(SIN(GAMA) ) ) 

LLEN  (L)  -  (CUH/S  IN  (THETA) )  -L2ET 

BLEN(L) -TPL-CUTL+CCUH/TAN ( THETA) )+TLOAD-(X*lE-01) 

LGAMA-SIN(THETA)*(XDELT/SIN(GAMA)) 

RLEN (L) -  (CLJH/COS(  2ETA) ) -LGAMA+JCL 
END  IF 

TIfC-Tire+1.0 
70  CONTINUE 

DO  75  M-l, LIMIT, 1 

PATH  (M)  -  CTT.EN  CM)  +LLEN  <M)  +RLEN (M)  +BLEN  (M) )  *lE-02 
75  CONTINUE 

K 

*  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  CHANGING  MAGNETIC  FIELDS. 

# 

IF  (TP7L.EQ.0.0)  THEN 
DO  80  N«l, LIMIT, l 

MAGH(N) "(TIMEI (N)*lE03)/(TPPU*lE-02) 

MAGB(N) -MAGH(N) *(4.0*3. !415927*lE-09) 

30  CONTINUE 

ELSE 

DO  85  P-l, LIMIT,  1 

MAGH(P) «(TIMEI  (P)*IE03)/LKP) 

MAGB(P)«MAGH(P) *(4.0*3. !415927*lE-09) 

35  CONTINUE 

END  IF 

$ 

*  THIS  PORTION  OF  THE  PROGRAM  USES  THE  STRENGTH  OF  THE  IDEAL  MAGNETIC 

*  FIELD  TO  CALCULATE  A  CHANGING  RESISTANCE  BASED  ON  THE  MAGNETIC  SKIN 

*  DEPTH.  IT  FURTHER  CALCULATES  A  RESISTANCE  DUE  TO  THE  RESISTANCE  OF  THE 

*  COPPER  PLATES.  IT  THEN  CALCULATES  THE  CHANGING  CURRENT  AND  MAGNETIC 

*  FIELDS  GENERATED  FROM  THIS  TIME  VARYING  RESISTANCE. 

* 
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*  MU  MAY  HAVE  TO  BE  CHANGED  IF  THE  METAL  IS  NOT  COPPER.  ROE0  DILL  HAVE 

*  TO  BE  CHANGED  IF  THE  METAL  IS  NOT  COPPER.  LITK  MAY  HAVE  TO  BE  CHANGED 

*  IF  THE  METAL  IS  NOT  COPPER.  CONS  LULL  HAVE  TO  BE  CHANGED  IF  THE  METAL 

*  IS  NOT  COPPER. 


a  £  a  X  IE  .1  .L  5,  a  a  A  A  A  a  a  a  a  a  .E  T" 


MU  "4. 0*3 . 1415927*lE-07 


LITK-1E-03 
ROE0-1.72*lE-08 
CONS- 11? 0.0 
DO  90  Q-l, LIMIT, 1 

TEMP(Q) »CONS#(MAGB(Q)#MAGB(Q)) 

ROE  (G) «ROE0*(  1 . 0+(L ITK*TEMP  CQ) ) ) 

SIGMft(Q) -1. 0/ROE (Q) 

KO (Q) -1 .0/(SIGMA(Q) *MU) 

CONTINUE 
THICK-TPT*lE-02 
DO  95  R=2, LIMIT. 1 
S-R-l 

TAU(R)  -MAGH(R)/(  (MAGH(R) -MAGH(S)) )/(  (MISEC(R)-MISEC(S))*1E-0S) 
SPHI(R)-2.0*SQRT(<i<O(R)*R*lE-0S)) 

CONTINUE 

DO  100  T-2. LIMIT. 1 

IF  (TPTL.EQ.0.0)  THEN 

REST  1  -TAU  <  T)  /  ( S I  GMA  ( T)  *MJ ) 

REST2 -SORT ( REST 1 ) 

REST«SIGMA(T)*REST2 
IF  (SPHI(T).LE. THICK)  THEN 

RES  (T)  -  ( 1 . 0/REST)  +(  (  ROE  (T)  *PATH  (T) )  /( ( TPPU*SPH  I  (T)  *lE-02) ) ) 
ELSE 

RES(T)  ■(  l .  0/REST) +(  ( ROE  ( T)  *PATH  ( T) )  /  ( TPPU*TH  ICK*1  E-02 ) ) 
END  IF 
ELSE 

RESTl-TAU(T)/(SIGMA(T)*MU> 

REST2 -SORT ( REST 1) 

REST -S IGMft ( T) *REST2 
IF(SPHI(T).LE. THICK)  THEN 

RES(T)-(1.0/REST)+((ROE(T)*PATH(T))/(W(T)*SPHI(T») 

ELSE 

RES(T)  -  <  1 . 0/REST)  +( (ROE  (T)  *PATH(T) )  /(IJ(T)  *TH  ICK)  ) 

END  IF 
END  IF 
CONTINUE 
RES (1) -RES (2) 

DO  105  U-2. LIMIT, 1 
V-U-l 

RAV(U)  - (RES (U) +RES ( V) )  /2 . 0 

LAV(U)  -  ( (TIMEL  (U)+T IMEL  (V) )  /2. 0)*lE-09+THROL+LOADL 
EARG 1 (U) - (RAV(U) /LAV(U) ) *- 1 . 0*lE-06 

EARG2  <U) ■ IO*(  (TIMEL ( 1 ) *lE-09) +THROL+LOADL)  /( (TIMEL (U)  *lE-09>  + 
+THROL+LOADL) 

CONTINUE 

EARG 1  ( l ) - (RES ( 1 ) / ( ( TIMEL ( 1 ) * 1 E-09) +THRQL+LOADL ) ) *- 1 . 0* IE-06 
EARG2(  D-IO 
DO  110  1-2. LIMIT. 1 
J-I-l 

EARG3( 1)>0.0 

EARG3  ( I )  -EARG  1(1)  4EA&G3  ( J) 

CONTINUE 

DO  115  K-l, LIMIT, 1 

EXPO (K) -EXP (EflRG3(K)) 


115  CONTINUE 

DO  120  L-l, LIMIT.  1 

IRES(L)-EARG2(L)*EXPO(L)*lE-03 
120  CONTINUE 

IF  (TPTL.EQ.8.0)  THEN 
DO  125  N-l, LIMIT.  1 

MAGRH  CM)  -  (IRES  (N )  ✓( TPPW*  IE-02 ) )  *  1 E03 
MAGRBCN) -MAGRH (N) *(4.0*3. 141 j927*1E-B9) 

125  CONTINUE 

ELSE 

DO  130  P-l, LIMIT. 1 

MAGRH  CP) - ( !RES(P)/U(P) ) *1E03 

MAGRB (P) -MAGRH (P) *<4. 8*3 . 14l5927*lt-09) 

130  CONTINUE 

END  IF 

* 

*  THIS  PORTION  OF  THE  PROGRAM  USES  THE  EMPIRICALLY  lEVELOPED  MODELS 

*  TO  DESCRIBE  THE  CURRENT  AND  MAGNETIC  FIELDS  WHICH  ACCOUNT  FOR  THE 

*  THE  FLUX  LOSSES  DUE  TO  A  SLIDING  CONTACT.  THE  LOSSES  DUE  TO  A  CHANGING 

*  AREA,  AND  THE  LOSSES  DUE  TO  A  FLUX  TRAPPING  EFFICIENCY. 

* 

OPEN  (UNIT-1. NAME-PLOT4,TYPE-'NEW') 

WRITE  (1.700) 

700  FORMRT  ('  PLOTS  OF  RESISTANCE  MODEL  CURRENTS  ') 

WRITE  (1,*)  LIMIT  . 

DO  800  I-I, LIMIT.  1 
J-I-l 

WRITE  (1.*)  J 
800  CONTINUE 

DO  805  K»l, LIMIT. 1 
A-IRES(K)*lE-03 
WRITE  (1.*)  A 
805  CONTINUE 

CLfcSE  (LiNIT-1) 

OPJN  (UN IT- l.NAft -PLOTS. TYPE-' NEW') 

URii'E  (1.705) 

705  FORMAT  ('  PLOTS  OF  IDEAL  MODEL  CURRENTS  *) 

WRITE  Cl.*)  LIMIT 
DO  810  L-l. LIMIT. I 
M-L-l 

WRITE  (1.*)  M 
810  CONTINUE 

DO  815  N-l. LIMIT. 1 
A-TIfEI (N)*lE-03 
WRITE  (1,*)  A 
815  CONTINUE 
* 

♦  THIS  PORI  ION  OF  THE  PROGRAM  SETS  THE  VALUES  FOR  THE  BEST  CASES  FOR 

♦  THE  EMPIRICAL  MODELS.  THIS  PORTION  OF  THE  PROGRAM  INITIALIZES  THE 

♦  SMALL  S  VALUES  FOR  THE  BEST  CASE  FOR  THE  SLIPPING  CONTACT  LOSS  MODEL. 

♦ 

II-LIMIT-18 
JJ -LIMIT-81 
AA-0.3102 
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BB *0.2015 
DEC1*BB/'(II-JJ) 

DEC2*(AA-BB)/'(LIMIT-II) 

* 

*  ■mis  PORTION  OF  THE  PROGRAM  INITIALIZES  THE  VALUES  FOR  THE  BEST  CASE 

*  EFFICIENCY  MODEL. 

* 

CKl-l 

KK-LIMIT-10 

LL-LIMIT-101 

CC-63.039 

DD-75.056 

DEC3- (DD-CC) /(L IMIT-KK) 

DEC4* ( 100 . 0-DD) /(KK-LL) 

* 

*  THIS  PORTION  OF  THE  PROGRAM  INITIALIZES  THE  VALUES  FOR  THE  BEST  CASE 

*  CHANGING  AREA  EFFICIENCY  MODEL. 

*  CO.. 
lii*LIMIT-5 
IJ-LIMIT-44 
IK-LIMIT-81 
EE -0.0009924 
FF *0.0006931 
GG *0.0005749 

DEC5- (EE-FF) x(L IMIT-MM) 

DECS* (FF-GG) /(Ml- I J) 

DEC7-GG/'(  IJ-IK) 

* 

*  THE  FIRST  ITERATION  PRODUCES  THE  BEST  CASE  PLOTS  FOR  THE  EMPIRICAL 

*  MODELS.  THE  SECOND  ITERATION  PRODUCES  THE  WORST  CASE  PLOTS  FOR  THE 

*  EMPIRICAL  MODELS. 

* 

OPEN  (UNIT* 1 . NAME-F ILNAM, TYPE-' NEW' ) 

WRITE  (1,710) 

710  FORMAT  ('  THE  FOLLOWING  RESULTS  AND  PLOTS  ARE  FOR  THE  BEST  CASE  ') 
IQ*IO*lE-03 
DO  900  1-1,2. 1 

* 

*  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  T^F  SMALL  S  VALUES. 

* 

DO  905  J-l.JJ.l 
SMS(J) -0.0 
CONTINUE 
JJ-JJ+1 
D1-DEC1 

DO  910  K-JJ.II.l 
SMS  OO  -D1 
D1-D1+DEC1 
CONTINUE 
D2-SM5(  1 1)  +DEC2 
U-II+I 

DO  915  L-I I.LIMIT, 1 
SMS(L) -D2 


905 


910 
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D2-D2+DEC2 
915  CONTINUE 

* 

K  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  EFFICIENCY  VALUES. 

DO  928  M-l.LL.l 
EFF(M) -100.0 
920  CONTINUE 

D3-DEC4 
LL-LL+1 

DO  925  N-LL.KK.l 
EFF(N) *100.0-D3 
D3-D3+DEC4 
925  CONTINUE 

IF  (CK1.EQ.1)  THEN 

D4- 100 . 0-EFF (KK) -DEC3 
KK-KK+1 

DO  930  P»KK, LIMIT, 1 
EFF(P)»180.8-D4 
D4-D4+DEC3 
930  CONTINUE 

ELSE 
END  IF 

DO  935  Q»l. LIMIT. 1 
EFF(Q)  •EFF(Q)/'100.0 
935  CONTINUE 
* 

*  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  CHANGING  AREA  FOR  THE 

*  CHANGING  AREA  EMPIRICAL  MODEL. 

* 

DO  948  R-l.IK,! 

AAA(R) "0.0 
948  CONTINUE 

IK-lK+1 
D5-DEC7 

DO  945  S-lK.IJ.l 
AAA(S) *D5 
D5*D5+DEC7  • 

945  CONTINUE 

IF  (CK2.EQ. 1)  THEN 
D6aAAA( IJ)+DEC6 
IJ-IJ+1 

DO  950  T-IJ.MM.l 
AAA(T) -D6 
D6-D6+DEC6 

958  CONTINUE 

D7*AAA  <MM) +DEC5 

itmth-i 

DO  955  U«MM,LIMIT, 1 
AAA(U) -D7 
D7-D7+DEC5 

955  CONTINUE 

ELSE 
END  IF 
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* 

THIS  PORTION  OF  TVIE  PROGRAM  CALCULATES  THE  CURRFNTS  DUE  TO  SUPPING 

*  CONTACT  LOSSES. 

$ 

DO  960  V«*l. LIMIT,  1 

IF (TPTL.EQ.0.0)  THEN 

POUER- 1 .0- ( (SMS (V) «2 . 0, /CUH) 

ISL  IP  < V)  -  IO*(  ( IRES  ( V)  /IO) )  *>*>OUER 
HSL IP  < V)  -  <  ISL  IP  (V)  /(TPPI»*lE-02) )  *IE03 
asup(V)  -HSLip(v>mj 
ELSE 

POUER  <*1 . 0-  ( (SMS  ( V)  *2 . 0)  /(U(  V)  *tE02) ) 

ISL  IP  (V)  -  1Q*(  ( IRES  ( V)  /IO)  )  **POUER 
HSL  IP  (V) «( ISL  IP  (V)  /U(  V) ) * IE03 
SSL  IP (V) «HSLIP(V)*MU 
END  IF 

'  060  CONTINUE 

•  * 

'  *  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  CURRENTS  DUE  TO  FLUX  TRAPPING 

*  EFFICIENCY. 

* 

DO  965  J-l. LIMIT.! 

IF  (TPTL.EQ.0.0)  THEN 
IEFF(J)-EFF(J)*IRES(J) 

HEFF(J)«(IEFF(J)/(TPPU*IE-02))*1E03 

BEFF(J)-HEFF(J)*MU 

ELSE 

IEFF(J)-EFF(J)*IRES(J) 

HEFF (J) • ( IEFF ( J )  Alt J) ) *!E03 
BEFF(J)-HEFF(J)*IU 
END  IF 
965  CONTINUE 
* 

*  THIS  PORTION  OF  THE  PROGRAM  CALCULATES  THE  CURRENTS  DUE  TO  THE  CHANGING 

*  AREA  EMPIRICAL  MODEL. 

* 

DO  970  KM. LIMIT.  I 
•  IF  (TPTL.EQ.0.0)  THEN 

LEX»( ( 1 .256637 l *AAA(K)  )/(TPPU*lE-02> )*lE-06 

A-EWOOO 

8-CNUM 

C  ■  ( (  T I MEL  ( K ) +L0 ADL ) * 1 E-09 )  +THROL-H.EX 
IAA(K)«IO*(B^C)*A 
HAA(K)«( IAA(K)/(TPPU*1E-02))*1E03 
BAA(K)-HAA(K)«MU 
ELSE 

LEX-( ( 1 . 25663?1*AAA (K) )/VJ(K) ) *lE-06 
A-EM>0(K> 

B-CNUM 

C  ■  ( (  T I  ft  L (K ) +L0 ADL ) * I E-09 ) +THROL+LEX 
IAA(K)«IO*(B/C)*A 
HAA(K)»(  IAAdOAKK)  )*!E03 
BAA(K) »HAA(K)*MU 
END  IF 
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970  CONTINUE 

URITE  <  1.7155 

715  FORMAT  ('  TIhC  IRES  ISLIP  IEFF  IAA#> 

DO  975  L-l, LIMIT, i 
M-L-l 

URITE  (1,720)  M» IRES(L), ISLIP(L) , IEFF(L) » IAA(L) 

720  FORMAT  (I6.G15.4,G15.4,G15.4,G15. 4) 

975  CONTINUE 

« 

*  THIS  PORTION  OF  THE  PROGRAM  SETS  UP  THE  PLOTTING  FILES  FOR  THE  SLIPPING 

*  CURRENT. 

* 

OPEN  (UNIT-1, NAPE -PLOTl. TYPE-* NEU') 

URITE  (1.580)  PLOT! 

URITE  Cl.*)  LIMIT 
DU  980  N-l. LIMIT. 1 
P-N-l 

URITE  (1.*)  P 
980  CONTINUE 

DO  985  Q-l, LIMIT, 1 
A- ISL IP(Q)  *1E~03 
URITE  (1,*)  A 
985  CONTINUE 

CLOSE  (UNIT-1) 

* 

*  THIS  PORTION  OF  THE  PROGRAM  SETS  UP  THE  PLOTTING  FILES  FOR  THE  EFFICIENCY 

*  CURRENT. 

* 

OPEN  (UNIT-1. NAIC-PLOTZ.TVPE-'NEU') 

URITE  (1,580)  PLOT2 
UR’TE  (1.*)  LIMIT 
D;>  990  R-l, LIMIT.  1 
S-R-l 

URITE  (1.*)  S 
998  CONTINUE 

DO  995  T-l. LIMIT, 1 
A-IEFF(T)*lE-03 
URITE  (1.*)  A 
995  CONTINUE 

CLOSE  (UNIT-1) 

* 

*  THIS  PORTION  OF  THE  PROGRAM  SETS  UP  THE  PLOTTING  FILES  FOR  THE  CHANGING 

*  AREA  CJRRENT. 

* 

OPEN  (UN  IT- 1 . NAIC-PL0T3,  TYPE-*  NELT  ) 

URITE  (1,500)  PLOT3 
URITE  (1.*)  LIMIT 
DO  1000  U-l, LIMIT, 1 
V-U-l 

URITE  Cl.*)  V 
1000  CONTINUE 

DO  1005  J-l. LIMIT. 1 
A-IAA(J)*lE-03 
URITE  (1,*)  A 
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1605  CONTINUE 

CLOSE  (UNIT-1) 

*  THIS  PORTION  OF  THE  PROGRAM  INITIALIZES  THE  VALUES  FOR  THE  UQRST  CftSE 

*  FOR  THE  EI’FIRICAL  MODELS.  THIS  PORTION  OF  THE  PROGRAM  INITIALIZES  THE 

*  S  VALUES  FOR  THE  SLIPPING  LOSS  MODEL. 

* 

I I-L  IMIT-45 
JJ-LIMIT-04 
AA-0.4800 
BB -0.1348 
DEC1-BBAII-JJ) 

DEC2-(AA-BB)/(LIMIT-II) 

*  THIS  PORTION  OF  THE  PROGRAM  INITIALIZES  THE  VALUES  FOR  THE  UORST  CASE 

*  FOR  THE  EFFICIENCY  MODEL. 

* 

CK1-0 
KK-LINIT 
LL-LIM IT-59 
CC -52. 8528 
DEC3-6.B 

DEC4-O80.0-CC)  ALIMIT-LL) 

DD-0.0 

*  THIS  PORTION  OF  THE  PROGRAM  INITIALIZES  THE  VALUES  FOR  THE  WORST  CASE 

*  FOR  THE  CHANGING  AREA  EtflRlCAL  MODEL. 

* 

CK2-0 

IJ-LIMIT 

IK-LIMIT-56 

EE-0.B0 18029 

DEC5-0.0 

DEC6-8.0 

DEC7-EEALIMIT-IK) 

FF-0.0 

GG-0.0 

OPEN  (UNlT-i.NArE-FNArC.TYPE-'NEVT) 

725  FORMAT* C:?2THE  FOLLOWING  ARE  THE  RESULTS  FOR  THE  UORST  CASE  ') 
9B0  CONTINUE 

i  THIS  PORTION  OF  THE  PROGRAM  PRINTS  THE  DATA  DEVELOPED  FOR  EACH 

#  OF  THE  MODELS. 

*  FORMAT  ('  THE  LINE  INDUCTANCE  IS:  '.F12.4,'  NANOHENRIES') 

P0rSaT(('2THE  GENERATOR  INDUCTANCE  IS:  '.F12.4.'  NANOHENRIES') 

FORMAT* ('2-THE  LOAD*"  INDUCTANCE  IS:  '.r 12.4.'  NANOHENRIES') 

URITE  (1.215)  CAP 

FORMAT<*  THE  BANK  CAPACITANCE  IS:  '.F12.4,'  MICROFARADS') 

URITE  (1,220)  VOLT 


200 

205 

210 

215 
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220  F0RP1AT('  THE  BANK  VOLTAGE  IS:  '.F12.4,'  KILOVOLTS' ! 

WRITE  (1,225)  EXPL 

225  FORMATS  THERE  ARE:  '.F12.4.'  KILOGRAMS  OF  EXPLOSIVES') 

WRITE  (1,230)  TAMW 

230  FORMATC  THERE  ARE:',F12.4,'  KILOGRAMS  OF  TAMPER' ) 

WRITE  (1,235)  MTOC 

235  FORMATS  THE  MASS  TO  EXPLOSIVE  RATIO  IS:  ',F12.4> 

WRITE  (1,240)  NTOC 

240  FORMAT('  THE  TAPPER  RATIO  IS:  '.F12.4) 

WRITE  (1,245)  0 

245  FORMATS  THE  DETONATION  VELOCITY  IS:  '.F12.4,'  PtVMICROSECONDS') 
WRITE  (1,250)  VA 

253  FORMAT('  THE  METAL  VELOCITY  IS:  ',F12.4,'  PtVMICROSECONDS') 

WRITE  (1,255)  THETD 

25  FORMAT('  THE  GURNEY  ANGLE  IS:  ',F12.4.'  DEGREES') 

WRITE  (1,260)  LIMIT 
260  FORMAT('  THE  LIMIT  IS:  ',110) 

WRITE  (1,265) 

265  FCRMAT('  THE  RESULTS  FOR  THE  FLUX  CONSERVED  MODEL  ARE:') 

WRITE  (1,270) 

270  FORMATS  TIME (MICROSECONDS)  INDUCTANCE (NANOHENRIES) * ) 

DO  150  T-l, LIMIT, l 

WRITE  (1,275)  MISEC(T) ,TIMEL(T) 

275  FORMAT('  ' ,4X,F12.4, 10X.F12.4) 

150  CONTINUE 

WRITE  (1,288) 

280  FORMAT ('  TIPE (MICROSECONDS)  CURRENT (K AMPS) ' ) 

DO  155  U-l, LIMIT. 1 

WRITE  (1,205)  MISEC(U).TIMEKU) 

285  FORMAT!'  '.4X.F12.4.6X.F12.4) 

155  CONTINUE 

WRITE  (1,290) 

290  FORMAT(*  TIPE  (MICROSECONDS) ' , 10X, 'HFIELD(APPS/METER) ' ,3X, 'BFIELD( 
+MEGAGAUSS) ' ) 

DO  168  V«l, LIMIT, 1 

WRITE  (1,295)  MISEC (V) , MAGH ( V) , MAGB ( V) 

295  FORMAT!'  ' ,4X.F12.4,4X,G28.4.2X,G20.4) 

160  CONTINUE 

WRITE  (1,300) 

300  FORMAT  ('  A  TIPE  VARYING  RESISTANCE  IS:') 

WRITE  (1,305) 

305  FORMAT  ('  TIME (MICROSECONDS )  RES I STANCE (OHMS) ') 

DO  165  1-1, LIMIT, 1 

WRITE  (1.310)  MISEC( I) ,RES( I) 

310  FORMAT('  ' ,F12.4.6X. G20.4) 

165  CONTINUE 

WRITE  (1,315) 

315  FORMAT  ('  THE  RESULTS  FOR  THE  TIME  VARYING  RESISTANCE  MODEL  ARE:') 
WRITE:  (1.320) 

320  FORMAT  ('  TIME (MICROSECONDS) ' ,3X, 'CURRENT(KAMPS) ' ) 

DO  170  J«l. LIMIT, 1 

WRITE  (1,325)  MISEC(J) , IRES(J) 

325  FORMAT  ('  ' ,4X,F12.4,6X,F12.4) 

170  CONTINUE 


83 


URITE  (1,330) 

330  FORMAT  ('  TINE(MICROSECQNDS)  *.  1 8X,  ’HFIELD(  AMPS/METER) ' ,3X. 'BFIELD( 
■HEGAGAUSS)*) 

DO  175  K-1, LIMIT, l 

URITE  (1,335)  MISEC (K) , MAGRHOO , MAGR8 (K) 

335  FORMAT  ('  ' ,4X.F12.4.4X,G28.4.2X.G20.4> 

175  CONTINUE 

URITE  (1,340) 

340  FORMAT  (*  THE  RESULTS  FOR  THE  SLIPPING  LOSS  MODEL  ARE : * ) 

URITE  (1,345) 

345  FORMAT  ('  T I  ME  ( M I CROSECONDS  )  * ,  3X*  *  CURRENT  ( K  AMPS )  * ) 

DO  180  L«l, LIMIT, 1 

URITE  (1.350)  MISEC(L) » ISLIP(L) 

358  FORMAT  ('  * .4X.F12.4.6X.F12.4) 

160  CONTINUE 

URITE  (1,355) 

355  FORMAT  ('  TIME(MICROSECONDS)M0X.'HFIELD(AMPS/METER)'.3X.'BFIELD( 
■HEGAGAUSS) 4 ) 

DO  185  M-l. LIMIT, 1 

URITE  (1.360)  MISEC (M) ,HSL IP (M) .BSLIP(M) 

360  FORMAT  ('  * ,4X.F12.4.4X,G20.4,2X,G2t.4) 

185  CONTINUE 

URITE  (1.365) 

365  FORMAT  (*  THE  RESULTS  FOR  THE  EFFICIENCY  MODEL  ARE:*) 

URITE  (1.378) 

370  FORMAT  ('  TIME (MICROSECONDS) * , 3X, *  CURRENT  (KAtES) * ) 

DO  198  N-l, LIMIT, 1 

URITE  (1.375)  MISEC (N).IEFF(N) 

375  FORMAT  (*  * .4X.F12.4.6X.F12.4) 

198  CONTINUE 

URITE  (1,380) 

380  FORMAT  ('  TIME  (MICROSECONDS)  * .  10X.  'HFIELDCAIES/TETER)  *  ,3X,  *BFIELD( 
■HEGAGAUSS)') 

DO  195  P-1. LIMIT. 1 

URITE  (1.385)  MISEC (P) ,HEFF(P) ,8EFF(P) 

385  FORMAT  ('  ' ,4X,F12.4,4X.G28.4,2X,G28.4> 

195  CONTINUE 

URITE  (1,390) 

390  FORMATS  THE  RESULTS  FOR  THE  CHANGING  AREN  LOSS  MODEL  ARE:') 

URITE  (1.395) 

395  FORMAT  (*  TINE  (MI CROSECONDS) ' ,3X, 'CURRENT (KAMPS)') 

DO  600  Q-l, LIMIT, 1 

URITE  (1,488)  MISEC (Q) , IAA(Q) 

400  FORMAT  ('  ' .4X.F12.4.6X.F12.4) 

600  CONTINUE 

URITE  (1,405) 

405  FORMAT  ('  TlfE(MICROSECONDS)',  l8X,'HFlELD(Aff»S/tETER)'.3X.'BFIELD( 
-HEGAGAUSS) ' ) 

DO  605  R-l, LIMIT, 1 

URITE  (1,410)  MISEC(R) ,HAA(R) ,BAA(R) 

410  FORMAT  (*  '  .4X.Fl2.4.4X.G20.4.2X.G2?.4) 

605  CONTINUE 
END 


84 


Appendix  B 


Appendix  B  conteine  the  curve*  of  the  aeasured  end  calculated 
data  uaed  for  analysis  in  Chapter  VI.  The  curve*  for  Teat  1  were 
presented  in  Chapter  VI.  This  appendix  contains  the  sane  curves  for 
the  regaining  tests.  They  include  the  experinental  currents,  the 
bubble  heights,  efficiencies,  and  additional  areas.  Also  included  are 
the  output  currents  of  each  of  these  models. 
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Figure  B-5.  Equivalent  Circuit  Current*  Test  2 


Figure  B-6.  Equivalent  Circuit  Currents  Test  3 
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Figure  B-8.  Equivalent  Circuit  Current*  Teat  5 
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Figure  B-10.  Time  Varying  Resistance  Teat 


Figure  B-12.  Tiae  Varying  Reaiarance  Test  5 
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Figure  B-15.  Slipping  Currents  Test 
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Figure  B-26*  Efficiency  Currents  Test  3 
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